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During the preheating era after inflation, resonant amplification of quantum field fluctuations takes 
place. Recently it has become clear that this must be accompanied by resonant amplification of 
scalar metric fluctuations, since the two are united by Einstein's equations. Furthermore, this 
"metric preheating" enhances particle production, and leads to gravitational rescattering effects 
even at linear order. In multi-field models with strong preheating (q 3> 1), metric perturbations 
are driven nonlinear, with the strongest amplification typically on super-Hubble scales (k — > 0). 
This amplification is causal, being due to the super-Hubble coherence of the inflaton condensate, 
and is accompanied by resonant growth of entropy perturbations. The amplification invalidates 
the use of the linearized Einstein field equations, irrespective of the amount of fine-tuning of the 
initial conditions. This has serious implications on all scales - from large-angle cosmic microwave 
background (CMB) anisotropies to primordial black holes. We investigate the (q, k) parameter space 
in a two-field model, and introduce the time to nonlinearity, t n i, as the timescale for the breakdown 
of the linearized Einstein equations. t n \ is a robust indicator of resonance behavior, showing the 
fine structure in q and k that one expects from a quasi-Floquet system, and we argue that t n i is a 
suitable generalization of the static Floquet index in an expanding universe. Backreaction effects are 
expected to shut down the linear resonances, but cannot remove the existing amplification, which 
threatens the viability of strong preheating when confronted with the CMB. Mode-mode coupling 
and turbulence tend to re-establish scale-invariance, but this process is limited by causality and for 
small k the primordial scale invariance of the spectrum may be destroyed. We discuss ways to escape 
the above conclusions, including secondary phases of inflation and preheating solely to fermions. The 
exclusion principle constrains the amplification of metric perturbations significantly. Finally we rank 
known classes of inflation from strongest (chaotic and strongly coupled hybrid inflation) to weakest 
(hidden sector, warm inflation), in terms of the distortion of the primordial spectrum due to these 
resonances in preheating. 
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I. INTRODUCTION 

The successes of the inflationary paradigm are based 
on classical perfection and quantum imperfections. In- 
flation leads, essentially because of the no-hair theorem, 
to a nearly maximally symmetric, de Sitter-like state. 
In contrast, it also leads, via quantum fluctuations, to 
small metric perturbations with a nearly scale-invariant 
spectrum Jj],||]. The subsequent evolution of these per- 
turbations, in the absence of entropy perturbations, is 
extremely simple — controlled by conserved quantities [^) . 
These allow one almost trivially to transfer the spectrum 
on super-Hubble scales from the end of inflation, typi- 
cally near the GUT scale ~ 10 16 GeV, down to photon 
decoupling where, in an f2 = 1 universe, the large-angle 
anisotropies in the cosmic microwave background (CMB) 
are formed, at a scale ~ 1 eV — which spans 25 orders of 
magnitude in energy. 

That such a very simple picture appears consistent 
with current CMB and large scale structure data is re- 
markable. In this paper we discuss some limitations of 
this picture and its fragility in the context of the recent 



revolution in inflationary reheating known as preheat- 
ing. We will show that for a wide range of models the 
post-preheating universe must necessarily have followed 
a much more complex route than the simple one outlined 
above. This more complex route typically invokes the 
full nonlinearity of the Einstein field equations and, in its 
most extreme guises, appears incompatible with current 
observational cosmology. In more moderate form, how- 
ever, the preheating route gives new freedom to the in- 
flationary paradigm in the form of its predictions for the 
CMB, large scale structure and primordial black holes. 

Recent study of the reheating era that follows inflation 
has highlighted a wide array of interesting new physical 
processes and phenomena. These new features all stem 
from nonperturbative, nonequilibrium effects in the infla- 
ton's decay. Whereas earlier treatments of post-inflation 
reheating modeled the inflaton as a collection of inde- 
pendent particles, each decaying into matter fields per- 
turbatively ||, these newer analyses focus on collective, 
coherent effects which may allow the inflaton condensate 
to transfer much of its energy rapidly into a highly non- 
thermal distribution of low- momentum particles. The in- 
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flaton's decay, in other words, may be resonant, coherent 
and dominated by the peculiarities of quantum statis- 
tics .0 During preheating, certain Fourier-modes of 
the matter-field fluctuations grow essentially exponen- 
tially with time, Sipk oc exp^^i), where jik is known 
as the Floquet index or characteristic (Lyapunov) expo- 
nent. Such rapid growth of the field fluctuations leads to 
particle production, with the number density of particles 
per mode k growing roughly as oc exp(2^£) [Q||U- 
The period of resonant inflaton decay has been dubbed 
"preheating," in order to distinguish it from earlier per- 
turbative studies of reheating. 

Details of this resonance, such as the rate of particle 
production and the momenta of particles produced, re- 
main model-dependent. Even so, certain consequences 
of such exponential growth in the occupation numbers of 
low-momentum quanta may be investigated. Such ef- 
fects include the possibility of non-thermal symmetry 
restoration and defect production supersymme- 
try breaking and GUT-scale baryogenesis jn}. These 
kinds of processes follow from the rapid amplification 
of matter-field fluctuations, driven during the preheat- 
ing era by the oscillation of the coherent inflaton con- 
densate. The strongly nonperturbative behavior within 
the matter-fields sector points to important cosmological 
consequences. 

In a recent paper Jl2| (hereafter Paper I), we high- 
lighted a distinct consequence which must in general arise 
during resonant post-inflation reheating: the resonant 
amplification of metric perturbations. Essentially, this 
arises because Einstein's equations directly couple the 
explosive growth in the field fluctuations to the fluctua- 
tions in the gravitational curvature, so that the previous 
approach of neglecting metric perturbations is inconsis- 
tent. Unlike some of the possible effects of preheating for 
the matter-fields sector, the rapid growth of metric per- 
turbations might leave distinct, observable signals in the 
CMB radiation, and could dramatically alter primordial 
black hole formation 

Here we develop the results of Paper I, via additional 
qualitative arguments and extensive numerical simula- 
tions of a two-field model. These simulations allow us 
to extend our discussion of the highly effective symbi- 
otic interaction between metric perturbations and par- 
ticle production during preheating. The source of the 
energy density for both metric perturbations and parti- 
cle production is the oscillating inflaton condensate. Be- 
cause the energy density of gravitational perturbations 
can carry an opposite sign to the energy density of the 
matter-field fluctuations, the two kinds of perturbations 
need not compete for the decaying inflaton energy, but 



* Loosely speaking, preheating is to the old theory of reheat- 
ing what the high-powered laser is to the household light-bulb. 



rather can grow resonantly together, coupled by the Ein- 
stein equations [Q. (Note that we consider only scalar 
modes of the metric perturbations here. An interesting 
extension would be to consider effects on vector or vor- 
ticity modes of the perturbed metric during preheating.) 

To emphasize this, contrast the resonant evolution of 
metric perturbations at preheating with that in tradi- 
tional studies of particle production in curved spacetime 
]L3[ . In the latter case, particle production arises via the 
polarization of the vacuum due to the shear of the met- 
ric and its perturbations, and hence energy conservation 
leads to a reduction in the anisotropy and inhomogene- 
ity of spacetimejj] Instead, the resonant amplification of 
metric and field fluctuations described here takes energy 
from the oscillating isotropic and homogeneous part of 
the gravitational field and so is qualitatively very differ- 
ent: both metric and field fluctuations can grow simulta- 
neously. 

In earlier studies of preheating, which neglected the 
presence and growth of metric perturbations, important 
distinctions were found between preheating in single-field 
models as opposed to models with distinct fields coupled 
to the inflaton . As highlighted in Paper I, and devel- 
oped further here, it is crucial to attend to these distinc- 
tions when studying the amplification of metric pertur- 
bations at preheating. When resonant decay is possible 
at all in the single-field cases, it is typically rather weak, 
occurring only within a narrow momentum range and 
governed by a small characteristic exponent fik- When 
distinct scalar fields are coupled to the inflaton, however, 
qualitatively different behavior results throughout most 
of parameter space: resonances may be very broad in mo- 
mentum space, and the rate of particle production may 
be much greater than in the single-field cases. Whereas 
the expansion of the universe during preheating is often 
sufficient to damp dramatically the resonant particle pro- 
duction in the single-field, narrow-resonance cases, broad 
resonances in the multi-field cases often survive in the ex- 
panding universe. 

As we demonstrate here, broad resonances in multi- 
field models can drive the exponential amplification of 
metric perturbations, and in these models, analyses 
based on Floquet theory remain useful, even when the 
expansion of the universe is included. (See Fig. (0).) 
In general, the effective characteristic exponents fik, c fi 
for the growth of metric fluctuations are greater than or 
equal to the exponents governing the field fluctuations' 
growth. 



' The most powerful use of this idea was in the Mis- 
ner "chaotic cosmology" program [jl4| whereby it was hoped 
that arbitrary initial anisotropy and inhomogeneity might be 
smoothed out by particle production JL5|. 
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FIG. 1. Growth of ln<i>fc versus mt, for q = 8000 with 
$ ~ 10~ 3 initially. $ is the gauge-invariant gravitational 
potential, m is the inflaton mass, and q is the resonance pa- 
rameter, as defined in Sec. |n| The specific model considered 
here is given by Eq. ([]). The main plot is for k = 0, and the 
inset shows growth for k = 20; the former is a super-Hubble 
mode, while the latter is within the Hubble radius. This plot 
comes from linearized perturbation theory, and so cannot be 
used to describe the modes' behavior in the nonlinear regime. 
Note, however, the strong breaking of scale-invariance in the 
earliest, linear-regime phase: the two modes go nonlinear at 
very different times. 

Furthermore, as is true of the nonthermal spectra of 
matter-field fluctuations amplified during preheating, the 
strongest resonances (largest fXk) for the metric pertur- 
bations occur for long- wavelength modes. Often these 
modes have wavelengths larger than the Hubble scale at 
preheating, k/a < H, where a(t) is the cosmic scale fac- 
tor. As discussed in Paper I, such long-wavelength modes 
may indeed be amplified causally at preheating; in gen- 
eral, this need not be surprising, since Einstein's field 
equations necessarily incorporate relativistic causality. 
The resonant growth of super-Hubble fluctuations derives 
from the coherence of the inflaton condensate on these 
large scales immediately after inflation, itself a natural 
consequence of inflation. This growth bypasses the stan- 
dard conservation law since it is accompanied by resonant 
growth of entropy perturbations, arising from the violent 
transfer of energy from the inflaton. Earlier studies of 
metric perturbations induced during (rather than after) 
inflation in multi-field models of inflation indicated that 
super-Hubble metric perturbations will not, in general, 
remain 'frozen' at their horizon-crossing amplitudes, but 
will rather evolve non-adiabatically p6| . The preheating, 
exponential amplification of super-Hubble metric pertur- 
bations studied here represents an extreme example of 
such dynamic, non-adiabatic, long-wavelength growth. 

By amplifying a distinct, highly scale-dependent spec- 
trum of metric perturbations after the end of inflation, 



multi-field preheating could dramatically alter the obser- 
vational consequences often assumed to hold for generic 
models of inflation. Rather than being sensitive only 
to the primordial, 'stretched' quantum fluctuations of 
the inflaton from deep within the Hubble radius during 
inflation, the power spectrum of density perturbations 
created during inflation would be deformed by the later 
preheating dynamics; the spectrum measured by COBE 
would therefore be sensitive only to these combined ef- 
fects. 

Given the robustness of the amplification of long- 
wavelength metric perturbations during strong preheat- 
ing, then, one of the greatest outstanding issues regarding 
the post-inflationary universe is the nonlinear evolution 
of these perturbations. The numerical simulations pre- 
sented here stem from integrating the equations of lin- 
earized perturbation theory. The time it takes for this 
linearized approximation to break down may be mea- 
sured self-consistently, and is discussed below. Typically, 
modes go nonlinear even before the matter-field fluctua- 
tions do, and hence well before what had previously been 
assumed to be the end of preheating. Beyond linear or- 
der, mode-mode coupling between metric perturbations 
of different wavelengths is likely to be crucial for under- 
standing the evolution of the power spectrum. These 
nonlinear effects would in general remove the individual 
'fingerprints' of specific inflationary models. 

Naturally, the challenge remains to study this fully 
nonlinear regime in quantitative detail. In preparation 
for such future studies, we indicate analytically below 
some effects which may be expected from such nonlinear- 
ities. We emphasize that these effects are qualitatively 
new. The old theory of perturbation evolution based on 
conserved quantities Q allowed for large amplification of 
metric perturbations at reheating, but required no pro- 
duction of entropy perturbations, i.e. purely adiabatic 
spectra on large scales. Perhaps more important, the con- 
served quantities were used to find the metric potential $ 
during inflation by using the COBE data AT/T ~ ICT 5 
and the Sachs- Wolfe relation AT/T w In the usual 
account, the large amplification of the k — * modes of $ 
between inflation and matter domination then requires 
fine-tuning of the inflaton mass (m ~ 10 _6 M p i) and self- 
coupling (A ~ 10 -12 ). 

What we discuss departs from this in two important 
ways. First, entropy perturbations are resonantly am- 
plified during preheating, so that the standard Bardeen 
parameter £ is not conserved and $ is further amplified. 
Second, and perhaps more important, we cannot use the 
total amplification naively to fine-tune the initial ampli- 
tude of <& during inflation to match with the $ ~ 10~ 5 
at decoupling required by the CMB. This second point is 
discussed further below, 
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in Sec. 

Work on the evolution of scalar metric perturbations 
through preheating began with study of the k — mode 
in simple single-field models P3,f8| , for which no super- 
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Hubble resonance was found. At this early stage, a major 
concern was the periodically singular nature of the evo- 
lution equations of the Bardeen potential $ due to the 
appearance of <j>~ x terms, where tp is the oscillating infla- 
ton background field. In contrast, the evolution of grav- 
itational waves (tensor perturbations) has none of these 
problems (l9|@. 

Historically, two-field models were studied soon after- 
wards 21 23]. In pH , the periodic singularities were re- 
moved by employing space-time averaging. This also re- 
moved the source of the parametric resonance and hence 
no parametric amplification was found. In p2| . the oscil- 
lations of the inflaton were retained and the singularities 
avoided by using the so-called Mukhanov variables. A 
resonance in the scalar metric perturbations was found. 
However that study was essentially restricted to classical 
mechanics - only the two background scalar fields were 
considered and no reheating to non-condensate modes 
(k > 0) was taken into account. Further, the study was 
limited to the first (and narrowest) resonance band of the 
Mathieu equation along the line A = 2q. Partly because 
of this, the authors never considered what happens when 
metric perturbations go nonlinear. 

Paper I |l^] was a response to the perceived limitations 
of these earlier works in regard to a realistic treatment of 
preheating. Work subsequent to Paper I has re-examined 
the single-field case |2^j24|] and begun study of the com- 
plex numerics of the problem when the metric perturba- 
tions become fully nonlinear ^5|. The need for such a 
program of study to understand strong preheating is one 
of the broad theses of this paper. 

The remainder of this paper is organized as follows. In 
Section |f| we review some basic features of preheating 
when the coupled metric perturbations are neglected al- 
together, and highlight the qualitative, as well as quanti- 
tative, differences between preheating with a single infla- 
ton field, versus preheating when the inflaton is coupled 
to distinct scalar fields. Section III presents a sketch 
of the important conceptual foundations of the alter- 
native picture of preheating introduced in Paper I and 
developed further here. Section |y] begins the detailed 
study of the coupled metric perturbations. We present 
the gauge-invariant coupled equations for fluctuations of 
the metric and fields, and discuss the issues of anisotropic 
stresses and entropy perturbations. In Section ^ we ad- 
dress the question of causality and the amplification of 
super-Hubble modes. 

In Section VI, we focus on a simple two-field model, 
performing a wide range of numerical simulations to 
show the growth and evolution of fluctuations, and the 
scale-dependence of amplification. The time t n \ during 
which the linearized approximation remains adequate for 
the resonantly-amplified metric perturbations is stud- 
ied in detail. Section VII turns to the essential role 



pling, and turbulence. We develop some of the conse- 
quences such strong nonlinearities may have for the per- 
turbation power spectrum, and discuss possible observa- 
tional consequences of the rapid amplification of long- 
wavelength metric perturbations. Section VIII discusses 



possible "escape routes" from the resonant amplification 
found in our simple model. In Section IX, we describe 



the hierarchy of inflationary models delineated by the 
strength of the amplification of super-Hubble metric per- 
turbations. Conclusions follow in Section |x| 



II. PREHEATING BASICS 
A. Minkowski spacetime 

In this section, we review briefly some of the key re- 
sults for preheating, as derived in the absence of metric 
perturbations. This section is meant to introduce key 
ideas and fix our notation; more complete reviews of pre- 
heating may be found in |J. We begin the discussion in 
Minkowski spacetime, and consider after that important 
deviations from these results for models in an expanding 
(but unperturbed) spacetime. 

The equation of motion for any scalar field \ coupled 
via a g4> 2 X 2 term to the inflaton, or for example, inflaton 
fluctuations in a model with a quartic self-coupling, takes 
the general form in Minkowski spacetime: 



£+nS(t) 



Xk = o, 



(i) 



Here . 



is the square of the time-independent 



natural frequency of the mode Xfc in the absence of in- 
teractions; ipo and uj are the amplitude and frequency 
of the inflaton's oscillations, respectively, and f(tut) = 
(j> 2 (t)/ipQ. The modes Xk thus evolve with a time- 
dependent frequency, flk(t)- 

The term £(i) measures backreaction on the modes Xk 
due to nonlinear x self-couplings, such as A x % 4 . A simi- 
lar term will appear in the equation of motion for the os- 
cillating inflaton field, which encodes the drain of energy 
from the oscillating inflaton to the amplified fluctuations. 
It is common in the preheating literature to treat such 
nonlinear terms with either the Hartree-Fock or large- N 
approximations, in w hich case £(i) ~ J d 3 k\xk\ 2 - As dis- 
cussed below in Sec. VII, these approximations capture 



played by nonlinear features in understanding the pre- 
heating epoch, such as backreaction, mode-mode cou- 



some elements of the nonlinear structure of the quantum 
theory, but neglect such effects as rescatterings between 
Xk modes with different momenta k. These rescattering 
terms can in many cases become quite important in the 
later stages of preheating. In any case, A x E(i) <C gtp% for 
early times, and this term may be neglected for the first 
stage of preheating. 
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Most models of inflation end with the inflaton <fi oscil- 
lating coherently. These include the common examples of 
chaotic inflation in convex potentials, and new inflation 
with its second-order phase transition, but excludes mod- 
els with first-order phase transitions which complete via 
bubble nucleation, warm inflation, and non-oscillating 
models. In those cases for which the inflaton does os- 
cillate coherently at the end of inflation, resonance ef- 
fects may become crucial. For early times after the start 
of the inflaton's oscillations, when A x £(t) is negligible, 
/(cut) is periodic when studying preheating in Minkowski 
spacetime. It is convenient to write 



n 2 k (t) =u? [A k +2qf(u;t)], 



(2) 



where Ak oc lu 2 /u> 2 , and q cx gip^/uj 2 . Exact relations 
depend on specifics of the inflaton potential V(4>); for a 
massive inflaton without self-couplings, q = gifig/ (Am 2 ), 
where m is the inflaton mass. (The factor of 4 in the 
denominator for q is a convention which comes when we 
write the equation of motion for Xk hi the form of the 
Mathieu equation.) In this early-time regime, the equa- 
tion of motion for the modes Xk is a linear second-order 
differential equation with periodic coefficients. Floquet's 
theorem then shows that solutions take the form |29l 



Xk(t) = Pk{t)e 



(3) 



where P k (t) is periodic. For certain values of fc, lying 
within discrete 'resonance bands' Afc, the Floquet index 
\x k will have a nonzero real part, indicating an exponen- 
tial instability. Modes within these resonance bands thus 
grow exponentially in time. For the remainder of this pa- 
per, we will write jj, k as the real part of the Floquet index 
only. 

Resonantly-amplified modes correspond to rapid par- 
ticle production. The occupation number-density for 
quanta of momentum fc may be written simply in the 
form 



&k 

2 



1 

m 



(4) 



This follows directly from standard Bogoliubov transfor- 
mation techniques |13f| ; its form may be understood in- 
tuitively, as discussed in ||, as = p k (x)/^k, where 
Pfe(x) = KlXfcl 2 + ^llXfcl 2 ) is the energy density carried 
by a given Xk mode. From Eq. (||), it is clear that modes 
within a resonance band, growing as Xk ~ exp(/Zfct), cor- 
respond to a number density of produced particles for 
that momentum of n k ~ exp(2/ifci). This is how preheat- 
ing produces a highly non-thermal spectrum of particles 
immediately after inflation: only particles in certain dis- 
crete bands Afc get produced resonantly, resulting in a 
highly non-Planckian distribution. Occupation numbers 
for quanta lying within the resonance bands grow quickly, 
and the small k modes may well be described as classical 



before the end of preheating, in the sense that the density 
matrix takes a diagonal form with rapidly oscillating off- 
diagonal terms which time-average (in Minkowski space- 
time) to zero and are suppressed in the long-wavelength 
limit d3JH]. 

The strength (/!&) and width (Afc) of the resonance 
depend sensitively upon the 'resonance parameter' q. In 
the narrow resonance regime, q < 1, analytic results yield 
fik ~ Afc ~ qui , indicating 'mild' exponential growth 
within narrow resonance bands ■ In the broad res- 

onance regime, q 1, analytic study of the coupled 
modes is far more involved, though significant progress 
has been made; in this case, the largest exponents sat- 
isfy (ik > 0(u>\, and the widths of the resonances scale as 
Afc oc q 1 I A y p|.p9|. In all known cases, the strongest res- 
onances (largest fj,k) occur on average for low- momentum 
(small-fc) modes. 

Still stronger resonances occur in models which pos- 
sess a 'negative coupling instability' [2S(] . This instability 
affects multi- field models for which the coupling g<t> 2 x 2 
between the inflaton and a distinct scalar field x is neg- 
ative: g < 0. (One can still maintain a positive-definite 
energy density for such models, and hence a stable vac- 
uum structure, by adding appropriate values of quartic 
self-couplings for both <fr and x-) Now the small-fc modes 
Xk 'feel' an inverted harmonic oscillator potential, that 
is, the parameters of Eq. (0) become approximately 



A k 



k 2 



in- 



2g, q 



9^1 
Am 2 ' 



(5) 



so that A h < for fc 2 



mi < 



Such models lead 



to highly efficient amplification of these modes Xk, now 
within even wider resonance bands than in the 'usual' 
broad-resonance cases. For modes subject to a negative 
coupling instability, the resonance band widths scale as 



Afc 



15 



1/2 



rather than as q 



1/4 



In Minkowski spacetime, the resonances which char- 
acterize the earliest stages of preheating will end when 
the backreaction of produced particles, E(t), significantly 
damps the inflaton's oscillations, ending the parametric 
resonance. This process is typically non-equilibrium and 
non-Mar kovian, i.e. it is very badly modeled by a T(f> 
term in the equation of motion for 4> |t],[| . One may find 
this time analytically, setting A x £(i cn d) = gifg and solv- 
ing, in the Hartree approximation, for t cna \ which yields 



to 



(4^n 



1 ]n(l0™mg 5 M pl 1 ), where /i max is the 



maximum value of fj, k . After t cna i the system of cou- 
pled oscillators enters the fully non-linear regime; even 
here, certain modes Xk may continue to grow due to 
strong rescattering effects ||| 30 . 



Before considering the important changes to these re- 
sults in an expanding, dynamical background spacetime, 
we first contrast preheating in a single-field model with 
preheating in multi-field models. Consider a single-field 
inflationary potential of the form 
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= W<? + ±A0 4 . 



(6) 



Considered in Minkowski spacetime, this model can pro- 
duce no inflaton quanta at all when A = 0: fluctuations 
5tp obey the same equation of motion as the background 
field </>, and simply oscillate. When A ^ and m ^ 0, 
the resonance parameter for the fluctuations Sip takes the 
form 



3A^| 
Am 2 



(7) 



Yet according to the usual treatment of the power spec- 
trum of primordial density perturbations generated dur- 
ing inflation |0,0|, both A and m for this model are 
constrained by observations of the CMB radiation to 
take the values A ~ 1CP 12 and m ~ 10 _6 M p i, where 
Mpi = 1.22 x 10 19 GeV is the Planck mass. Furthermore, 
for chaotic inflation, the slow-roll phase will end, and the 
inflaton will begin to oscillate, at ipo ~ 0.3M p i. These 
yield q < 0.1. When m — 0, the inflaton's frequency of 
oscillation is well-approximated by oj « 0.85-s/A^o [pLQlj 
yielding q ss l.[j] Needless to say, smaller values of the 
inflaton's amplitude, corresponding for example to GUT- 
scale inflation, will of course yield far smaller values of q. 
Thus the maximum values of the resonance parameter, 
q, in single-field models of inflation are restricted to the 
narrow-resonance regime, q < 1. 

The restriction to the narrow-resonance regime for 
single-field models of inflation may be contrasted with 
the situation in multi-field models. Consider, for exam- 
ple, the simple model 



V(<t>,x) = km 2 <j> 2 + h<f> 2 X 2 



(8) 



In this case, the resonance parameter for the x field be- 
comes 



q 



9Vo 
Am? 



gx 10 



10 



(9) 



where the last expression on the right holds for the 
chaotic inflation case. In this case, CMB observations 
only restrict the coupling g through radiative loop cor- 
rections, which suggest \g\ < 10 6 -^] Thus even weak 



couplings between the inflaton and the distinct scalar 
field x wu l produce resonance parameters well into the 
broad- resonance regime, with \q\ 1. Multi- field models 
of preheating generically lie within the broad-resonance 
regime, whereas single-field models are necessarily lim- 
ited to the narrow-resonance regime. This difference will 
be crucial when we consider the amplification of metric 
perturbations at preheating later. 

To finish the discussion of preheating in Minkowski 
spacetime, we examine the case of non-periodic inflaton 
evolution. Perhaps the simplest and most unified way to 
do this is in terms of the 1 — D Schrodinger operator, C, 
which is dual to the Klein-Gordon operator in Eq. (|l|) 
under the interchange of time and space p3] . Under this 
exchange, the unstable modes k for which fj,k > are 
identified with the complement of a subset of the spec- 
trum of £ (technically the complement of the absolutely 
continuous part of the spectrum, oac p2fl )- From this 
we immediately see why always had a band structure 
in the cases considered above, since in the dual picture 
we simply find the conduction bands of condensed matter 
models of metals with periodic potentials - the Brillouin 
zones. 

In the case where the inflaton evolution is quasi- 
periodic, the above duality enables one to show that the 
band structure of the periodic case is destroyed, leav- 
ing only a nowhere-dense, Cantor set of stable modes 
|52j,[3l]]. Similarly, in the case that 4> exhibits partial or 
complete stochastic behaviour, one finds that all modes 
grow (/ife > 0), with measure one p2| . In the dual, con- 
densed matter picture, this simply corresponds to the 
famous Anderson localization of electron eigenfunctions 
in metals with random disorder [ B3| . This was found in- 
dependently by Zanchin et al. p4|, who also showed that 
it leads to an increase in fj.^ over the purely periodic case 
at the same q, i.e. noise enhances preheating. 

An important application of these results is the case 
when (f> evolves chaotically in time, as is expected when 
there are many fields coupled to the inflaton J35[. In the 
strongly-coupled limit one can show that the <j) motion 
becomes stochastic and hence one knows that preheating 
is enhanced in this region of the parameter space |3l[| . 



"'When the inflaton potential contains quartic self-couplings, 
the specific resonance structure for coupled fields is quite dif- 
ferent from the case when A = 0; in particular, whereas an 
infinite hierarchy of resonance bands results when A — 0, only 
one single resonance band exists for A 7^ 0, m = 0. Even 
so, in both cases, preheating in a single-field model is restricted 
to the narrow-resonance regime, for which q < 1 and /i^ ~ qcu, 
and thus these fine-structure points are not relevant to our 
present discussion. 

s This constraint does not exist in (softly-broken) supersym- 
metric theories. Hence the inflaton may be strongly coupled 
to other fields and the resonance parameters may be arbitrar- 
Uy large El. 



B. Unperturbed FRW background 

Now consider what happens for preheating in an ex- 
panding FRW (unperturbed) universe with scale factor 
a(t), rather than in Minkowski spacetime. In this case, 
a generic equation of motion for modes coupled to the 
inflaton takes the form 



dt 2 ^ 



n 2 k (t) 

k 2 
= — ■+ 



0. 
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-H 2 



3 a 
2 a 



A x E(t), 



where H = a/a is the Hubble constant. Note that both 
the physical momentum, k/a(t), and the amplitude of the 
oscillating inflaton, </?o(i)> now redshift with the expan- 
sion of the universe. Note also that Floquet's theorem 
no longer applies, since the terms within 0?(t) are no 
longer simply periodic. This translates into the phys- 
ical result that no resonant amplification occurs in an 
expanding universe for models restricted to the narrow- 
resonance regime. In particular, single-field models of 
inflation which include a massive inflaton fail to produce 
any highly-efficient, resonant particle production at the 
end of inflation, when the expansion of the universe is 
included .Q 

Explosive particle production is still possible in an ex- 
panding universe, however, for multi-field models in the 
broad-resonance regime, with q > 10 3 |H . (For models 
in which the \ field has a mass comparable with the infla- 
ton's mass, resonant growth occurs for q > I0 4 .) In this 
case, the growth of the \k modes does not take a sim- 
ple exponential form; instead, certain modes Xk jump 
stochastically each time the oscillating inflaton passes 
through zero, but evolve quasi-adiabatically in between 
these sudden jumps ||. To understand these jumps, it is 
convenient to work in terms of an adiabaticity parameter, 



«2 



The large amplification of the modes Xk occurs during 
those periods when r > 1. For definitcness, consider 
again the model of Eq. (||). In an expanding universe, 
for early times after the start of the inflaton's oscillations, 
the inflaton will evolve as 

4>{t) = <p (t) eun(mt) « y 0) \ (12) 

mt 

and an effective resonance parameter for the Xk modes 
may be written 



q e ft{t) = 



Am 2 



(13) 



Then, for q e g 3> 1, we have, to leading order in t , 

-1/2 



g cff cos(mt) 
2 sin 2 (mi) 



(14) 



** An exception comes with the model of a massless, quarti- 
cally self-coupled inflaton which is conformally coupled to the 
curvature, for which weak resonance of the Sip fluctuations 
survives in an expanding universe, again with Uk_^ qto < lo, 
due to the conformal invariance of the theory 



(10) Given that q e g 3> 1, we see that we will only get strongly 
non-adiabatic evolution of the modes Xk near the zeros 
of 4>{t). Analytic estimates exist for fi). in this case [||. 
Two points are of special interest. First, /i^ is largest for 
k = 0, with /i fc oc ln(l + 2 cxp[-vrfc 2 /(2a 2 m 2 ^i)] + V). 
Hence the amplification of super-Hubble modes is present 
even in the absense of metric perturbations. Second, /xj, 
depends on the phase of the wavefunction of %, through 
the term V - a purely quantum effect. The change in the 
phase between the zeros of 4> is much larger than 2ir, and 
hence the phase (mod 2ir) is a pseudo-random variable 
of time. This stochastic behavior is transmitted to 

Thus, for the broad-resonance regime, strong ampli- 
fication is still possible at preheating, though it is no 
longer of a simple parametric resonance kind. Never- 
theless, many concepts from the pure Floquet theory re- 
main helpful, in particular q c g and (J>k,es(t) = Xk Xk 
fyipsll . As in the Minkowski spacetime case, for q c g ^S> 
1, the effective characteristic exponents may still be 
large, /ik,es(t) > lo, though now they are explicitly time- 
dependent. One way to understand the evolution of the 
amplified Xk modes in this case is that as a given mode, 
k/a(t), redshifts with the expanding universe, it will pass 
through a series of broad resonance bands, Afc, and get 
amplified in each one ^,^6|. As in Minkowski spacetime, 
the widths of these resonance bands for q D g 3> 1 scale as 
Afc ~ q^ 4 ||. Similarly, the widths of resonance bands 
for models with a strong negative-coupling instability in 

(11) an expanding universe scale as Afc ~ | ^eff | x / 2 [f29f . 

In the old theory of preheating, resonant passing of 
energy back to the inflaton only occurs at second order 
through rescattering effects. This was first discovered nu- 
merically H and then explained analytically || as arising 
from mode-mode coupling, which acts as a driving term 
that leads to population of higher-momentum (shorter- 
wavelength) modes of the inflaton. This provides a way 
of stemming the build-up of backreaction effects and al- 
lows the total variances of the fields to be larger than 
would be estimated using only the Hartree or N — > oo 
approximations, which miss all rescattering effects. This 
is of relevance to the discussion of the maximum scale for 
non-thermal symmetry restoration [B 10 1. In Sec. IV we 
show that when the metric perturbations are included, 
they induce rescattering even at linear order and again 
this helps to raise the maximum variances. 

Multi- field models, themselves inherently more real- 
istic than models which include only an inflaton field, 
lie generically within the broad-resonance regime, with 
q e g » 1. Resonant growth of these coupled field modes, 
then, marks the immediate post-inflation era, as preheat- 
ing produces large occupation numbers for quanta in 
various resonance bands Afc. The highly non-thermal 
distribution of particles results from a stochastic process, 
in which modes Xk slide through a series of resonance 
bands and receive a large, non-adiabatic 'kick' each time 
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the oscillating inflaton passes through a turning point. In 
the broad-resonance regime, this strongly non-adiabatic, 
rapidly-growing behavior of certain modes Xfc remains 
robust even when expansion of the universe is included. 

As the remaining discussion will demonstrate, it is 
not only field modes Xk which enjoy such robust reso- 
nant behavior at preheating; metric perturbations do as 
well. Moreover, the amplified metric perturbations in 
turn provide one more source to induce increased par- 
ticle production, by means of gravitational rescattering. 
By approaching the immediate post-inflation epoch in a 
fully self-consistent manner, the surprising behavior of 
preheating becomes all the more striking. 

III. CONCEPTUAL FOUNDATIONS OF THE 
METRIC PERTURBATION RESONANCES 

Before we embark on a detailed study of the ampli- 
fication of metric perturbations during preheating, it is 
perhaps appropriate to highlight the important concep- 
tual differences between metric perturbation evolution in 
the case of a single field and the case involving many 
fields interacting both via gravity and other forces.^ 

As we will see below (see Eq. (ffl|)), Einstein's field 
equations couple metric perturbations directly to field 
fluctuations Sipik, for all wavenumbers k. This means 
that metric perturbations can in general be amplified 
at any k. On small scales, k/a > H, such resonances 
could lead to nonlinear effects such as pattern and pri- 
mordial black hole formation [ |37|]38| ]. On much larger 
scales, k/a <C H, there is in addition the question of the 
causal amplification of such long- wavelength modes 

The schematic in Fig. (|^) shows the five main points 
which are involved in a conceptual understanding of met- 
ric perturbation resonances in strong preheating. All ex- 
cept causality are related to the existence of multiple 
fields in realistic preheating models, and are relevant for 
large-fc and small-fc resonances alike. 

• Non- gravitational forces may dominate 

In the case of fields interacting through non- 
gravitational forces, the coupling constants may 
naturally be large, implying that the evolution of 
metric perturbations is non-gravitationally domi- 
nated |3l]]. This can destroy the solution $ = 
const, for the k — growing mode. 

• Large resonance parameters 



T^The case of many non-interacting fields is well studied, but 
excludes preheating, and is typically restricted to slow-roll. 
See e.g. @. 



In the single-field case, the natural dimensionless 
resonance parameter q is order unity or smaller 
since both the coupling and frequency of oscilla- 
tions is controlled by the same constant of the the- 
ory (e.g. the effective mass or self-coupling). In the 
multi-field case, the couplings between fields may 
differ from the inflaton frequency of oscillation m 
by orders of magnitude, and hence the associated 
resonance parameters qi = gi^p 2 /m 2 can naturally 
be very large. 

• Entropy perturbations 

In the multi-field case, entropy perturbations are 
generic (see discussion above Eq. (^5|)). In the 
presence of strong non-gravitational interactions 
and the associated explosive transfer of energy be- 
tween fields with different effective sound speeds, 
these entropy perturbations are resonantly ampli- 
fied, thereby invalidating the use of standard con- 
served quantities to calculate <&. 

• Causality is not violated 

Causality is not violated by resonant amplification 
of the k ~ modes. This is discussed further in 
Sec. [y| in terms of the unequal-time correlation 
function for the perturbed energy density. The 
field equations are relativistic and hence causality 
is built into their solution. The apparent violation 
of causality arises essentially from the 'initial con- 
ditions' at the start of preheating: the homogene- 
ity and coherence of the zero-mode of the inflaton, 
which is required to solve the horizon problem. 

• Multi-fields 

/ Non-gravitational forces 

Super-Hubble I 

Resonances *T •Causality 

\ _ Large resonance parameters 

V 

Entropy perturbations 

FIG. 2. Conceptual points of importance in understand- 
ing the origin of the metric perturbation resonances in strong 
preheating. 

Note that here and throughout this paper, we use the 
term "super-Hubble" to refer to modes with wavelengths 
longer than the Hubble radius, that is, with k/a < H . It 
is common, though somewhat misleading, to use the term 
"super horizon" for this. As is well-known, the Hubble 
radius marks the speed-of-light sphere in any given cos- 
mological epoch, but can be dramatically different from 
the actual particle horizon. Of course, by solving the 
horizon problem, inflation gcncrically produces a par- 
ticle horizon which is many orders of magnitude larger 
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than the Hubble radius during post-inflationary periods 
in the universe's history. As emphasized in this section 
and below, the "super-Hubble" resonances described in 
this paper are thus amplified fully causally; their wave- 
lengths remain smaller than the true particle horizon af- 
ter inflation, though they stretch beyond the Hubble ra- 
dius at the end of inflation. Related points have been 
made in the study of the backreaction of gravitational 
waves H^]: not only are the super-Hubble modes con- 
sistent with causality and locality, but if one imposes a 
Hubble-scale infrared cut-off and excludes these modes, 
then energy conservation is violated. 

Throughout we simply write k = for the longest- 
wavelength modes which we consider. Physically, the 
maximum wavelength of perturbations which can ever 
be amplified during preheating corresponds to the cor- 
relation length of the coherent inflaton condensate. In- 
flation guarantees that the inflaton background field will 
be coherent (and thus will oscillate at preheating with a 
spacetime-independent phase) over vastly larger length 
scales than the Hubble radius at preheating. This is 
what explains the causal amplification of super-Hubble 
perturbations. Thus, where we write k = 0, we really 
intend something more like k ~ 10 -28 , corresponding to 
length-scales well within the Hubble radius during infla- 
tion, which then were stretched coherently during the in- 
flationary epoch. Of course, in our numerical simulations, 
these exponentially-tiny wavenumbers remain practically 
indistinguishable from 0. We return to the question of 
the inflaton condensate's coherence length below, in Sec. 



VII. 



IV. METRIC PREHEATING 

With the points of the previous section clearly in view, 
we may proceed to study in detail the behavior of met- 
ric perturbations during preheating. We use the gauge- 
invariant formalism of |2| to investigate the evolution of 
scalar perturbations of the metric. (The amplification of 
tensor perturbations, or gravitational waves, at preheat- 
ing has been considered in (l9| .) In longitudinal gauge, 
the perturbed metric takes a particularly simple form, 
expressed in terms of gauge-invariant quantities: 

ds 2 = (I + 2$) dt 2 - a 2 (t) (1 - 2*) 5 lJ dx i dx j , (15) 

where $ and \& are gauge- invariant potentials. On length 
scales shorter than the curvature scale (which itself is 
longer than the Hubble length), preheating in an open 
universe proceeds with little difference from the spatially- 
flat case (see the papers by Kaiser in 0), so we restrict 
attention in this paper to the case of a spatially-flat back- 
ground. 



A. The single-field case 

To begin, consider a single-field model, with the in- 
flaton field separated into a background field and quan- 
tum fluctuations, (f)(x) = ip(t) + 5(p(x), where Sep is the 
gauge-invariant fluctuation in the longitudinal gauge [0 . 
The Friedmann and Klein-Gordon equations governing 
the background quantities are 



H 2 = \n 2 + V(<p)] , 
<p> + 3Hcp + V v = , 

where V v = dV/dip and k 2 = 8ttG = SnM- 2 
may be combined to give: 

H 

The energy-momentum tensor 

T» v = WV"(/> - g" v [iV Q 0V a - V(<f>)] 
may be written in the form 



12-2 

-2 K V ' 



T" v = {p + p)u^u v 



pg 



(16) 
(17) 

These 



(18) 



(19) 



(20) 



where u M oc V^(/> is the natural four- velocity defined by 
the scalar field, which coincides in the background with 
the preferred 4- velocity. The energy density and effective 
pressure are 
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1 12 



+ V, p- 



V . 



where an overdot denotes the comoving covariant deriva- 
tive m^V^j. Thus in the natural frame, T^ v has perfect 
fluid form, so that the energy flux and anisotropic stress 



vanish in this frame: = 
which is close to u^, i.e., u M = 
we have to linear order 



= 7i), v . In any other frame 
u v _|_ v m where v u u^ = 0, 



P-, p = p, q^ = -(p + pH 



0. 



(21) 



Thus, to linear order, the scalar field has no anisotropic 
stress in any frame that is moving non-relativistically rel- 
ative to the natural frame. Writing T^ v = T^ v + 5T^ V , 
where T^„ is the background energy-momentum tensor 
and 5T^ V is the (gauge-invariant) perturbation, it follows 
that 



ST l j oc 8 % j 



(22) 



Then the perturbation of the Einstein equations SG^ = 
k 2 ST u 



fiu gives 



= for i^j. 



(23) 



so that * = $. 

The gauge-invariant perturbations for this single-field 
case then satisfy the equations 
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3#$ + 3H 2 - V 2 $ 



\k 2 



<p(6<p)' + V ip Stp- tp 2 § 



{Stp)' +3H(5tp) 
4<p$-2V v §, 
3' 



(24) 

(25) 
(26) 



Note that the metric perturbations enter at the same, 
linear order in Eq. ( |25| ) as the usual term, V w 8tp, and 
must be included in this equation to remain consistent. 
Eq. (E9) can be integrated directly to give: 



2a 



atpStp dt . 



(27) 



which allows us to eliminate $ from Eq. fl25|) at the 
expense of introducing terms which depend on the com- 
plete time-history of tp and dtp. In the realistic case where 
the fields exhibit a stochastic component to their evolu- 
tion 



HM32 34 1, this non-Markovian term in the evolution 
equation is expected to cause new conceptual deviations 
from the simple Markovian equation in the absence of $ 



By Eqs. ( |19[ ) and (|20|), the perturbed energy density 
and pressure in longitudinal gauge are 



Sp = tp(Sp)' + V<p6tp - tp 2 § . 
Sp = tp(5p)' — VipStp — tp 2 <& . 



(28) 
(29) 



Eq. ( p8| ) shows how the energy density carried by grav- 
itational fluctuations, 6p[<&\ = —tp 2 ^, is negative when 
$ > The entropy perturbation Y is defined by 



c 



2H 

K 2 tp 2 



(32) 



in the single-field case. It follows from Eqs. @, ©, 
© and (H) that 



c 



a 2 pp 



Provided that tp (and hence H) is always nonzero, the 
k = Q mode of ( is conserved for models with only one 
scalar field, and the k — entropy perturbation is zero. 
However, during reheating, tp periodically passes through 
zero, and this can produce entropy perturbations even on 
super-Hubble scales, as noted in jl8|,^3]]. These entropy 
perturbations will not be strong, since at most weak res- 
onance occurs on super-Hubble scales in single-field mod- 
els. 



B. The multi-field case 

We now generalize the above discussion and study 
the behavior of <& for models with N interacting scalar 
fields (j>i, where I = 1,---,N, taking the inflaton to 
be I = 1. Each field may be split into a homoge- 
neous part and gauge- invariant fluctuations as (f>i(t, x) = 
tpjit) + 6tpi(t,x). The background equations, Eqs. Jl6|), 
(Eg), become 



1 ~2 



H 



(33) 
(34) 



pT = 5p — — } 5p . 



Using Eqs. (E4)-(E9J), we can show that 



pT 



414 



3n 2 a 2 Htp 



V 2 $. 



(30) 



Thus on all scales k > 0, the single scalar field generi- 
cally induces nonzero entropy perturbations. These are 
easily missed if one neglects the metric perturbations 
Slow-roll inflationary conditions will keep these entropy 
perturbations negligible during inflation, but they will 
grow in concert with the metric perturbations $ during 
the oscillatory, preheating phase. 

An important quantity is the Bardeen parameter, i.e., 
the spatial curvature on uniform density surfaces M 



H 



(31) 



which evolves as 



The first question to address in the multi-field models is 
whether anisotropic stress can arise at linear level. Each 
field defines its natural frame via Uj oc V^j, and these 
are all close to each other in a perturbed FRW universe, 
and all reduce to the preferred frame in the background. 
Given a choice u M of global frame that is close to these 
frames, we have Uj = + Vj , where the relative veloc- 
ities satisfy VjU^ — 0. The energy-momentum tensor in 
Eq. ( |l9| ) becomes 



9" 



and Eq. (|2C|) generalizes to 

= {p + p)u^u" - pg^ + u»q v + g"u" , 

where the total energy density, pressure and energy flux 
are 
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IE- 
i 



E i2 ^ 
i 



and the anisotropic stress vanishes at linear order. 
Since these total quantities obey the same frame- 
transformation laws as in Eq. (|2l|), it follows that the 
anisotropic stress vanishes to first order in any frame 
close to u^. Then, as before, Eqs. ( f22| ) and ( p3| ) hold, 
and we have VP = <E>. We emphasize that the vanishing of 
anisotropic stresses in a multi-field universe holds only in 
linearized gravity (see |4lf for the nonlinear equations). 
If fluctuations grow to nonlinear levels and relative ve- 
locities become relativistic, then anisotropic stresses will 
emerge and could significantly affect the dynamics (com- 
parejH). 

The nature of entropy perturbations is another exam- 
ple where the single-field and multi-field cases bifurcate 
strongly. In the multi-field case, it is known that en- 
tropy perturbations are generic, arising from differences 
in the effective sound speeds, from relative velocities, and 
from interactions (see, e.g., p0| , p3[ ). Since reheating is 
the conversion of energy from a coherent field to a rela- 
tivistic, thermalized gas, it is clear that entropy will be 
generated. Even exactly adiabatic initial conditions at 
the exit from inflation will not prevent the production of 
entropy perturbations during reheating, since the fields 
are interacting non-gravitationally. 

The entropy perturbation in a multi-component sys- 
tem is defined by Effl 



P r = E ( 6pi ~ c s S Pi) ' 
i 

where the 'total' effective sound speed is given by 

< 2 4e^ 2 
i 

where hj = pj + pi , h = hi . 



(35) 



Then T can be related to (. 



C = -^V 2 $ + -H 
H 2 



In the two-field case, (, is 



r * 2 • 2 - 


( ^fl 




.01+02. 


V 01 


02 / 



which shows explicitly how Q can grow even for k — 0. 
As we have discussed in Sec. [n], preheating with more 
than one scalar field is a period of intensely non-adiabatic 
evolution, and we may expect that £ will deviate strongly 
from even in the k — > limit. This is precisely the kind 
of non-adiabatic growth which we observe below. This 
invalidates the use of £ as a way of finding $, and one is 
forced to explicitly integrate the evolution equations to 
find <&, as we do in Sec. VI. Further discussion of these 
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as we do in Sec 
points is made in Sec 

Note also that £ may be written in terms of the 
Mukhanov perturbative quantities 



Qi = 8<pi + 



In the two-field case [E2[l 



c 



H {iplQl + </32<?2) 



0i 



02 



(37) 



(38) 



The linearized equations of motion for the Fourier 
modes of 3> are (see where we have corrected the 
first equation) 

3H$ k + [(k 2 /a 2 )+3H 2 ] $ k = 

-5« 2 E [pi( 6c fn)' -^j + VjStpjk , (39) 
i 

(Sip Ik )" + 3H (Sip Ik )' + (k 2 /a 2 )Sip Ik = 

4$ fe 0j - 2V^ k - V ™ 5 VJh, (40) 
j 

$ fe + H$ k = \k 2 <Pio~<fik ■ (41) 

Equation ( [ill ) exhibits a simple generalization of the inte- 
gral solution of the single field case, Eq. (^7|). There are, 
however, important distinctions to appreciate between 
the single-field and multi-field cases. Note first that when 
more than one scalar field exists, cross-terms arise in the 
equations of motion for each of the field fluctuations, Sifj. 
That is, the evolution of Sipj depends, even at linear or- 
der, on all the other Sipj. These cross-terms have been 
neglected in previous studies of multi-field models of pre- 
heating, which ignored metric perturbations |^]-[|. They 
are present here, in general, because the gravitational po- 
tential $ couples universally with all matter fields, that 
is, with the entire ST^ U . Each matter field fluctuation 
6(pi 'talks' with the metric perturbations, which in 
turn 'talks' with all of the other matter field fluctua- 
tions. Given the form of Einstein's equations, then, each 
of these field fluctuations 5ipi becomes coupled to each 
of the others. 

Though distinct from nonlinear mode-mode coupling 
of the sort discussed in Sec. VI], these cross-terms, enter- 



ing at linear order, can play a large role in making "grav- 
itational rescattering" quite effective at preheating. We 
adopt the term rescattering for these linear-order cross- 
terms, whereas it is usually reserved for discussing non- 
linear effects, because the coupled form of the equations 
in Eqs. (^9|)-(|4l|) indicates that modes are amplified, 
and quanta are produced, due to the coupling between 
fluctuations of different fields. The significance of these 
linear-order rescattering terms is confirmed in our nu- 
merical simulations, presented below. 

As discussed in Sec. |l[ further differences remain 
generic between single-field and multi-field models at pre- 
heating, and these distinctions result in dramatic dif- 
ferences for the metric perturbations in each case. In 
the single-field case, resonances, when they survive the 
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expansion of the universe at all, are quite weak. Fur- 
thermore, the amplitude of the inflaton's oscillations de- 
cays with monotonic envelope due both to expansion of 
the universe and to dissipation from particle production. 
The right hand side of the constraint equation, Eq. (26), 



end this is not crucial for the robustness of the resonances 



is therefore either weakly growing or decaying in time. 
Hence in this single-field case is either weakly grow- 
ing or decaying in time. 

In the multi-field cases, we may expect far greater res- 
onance parameters, q c g 3> 1, resulting in much larger 
Mfc,cff than in the single-field case. That is, in the multi- 
field case, the system will in general be dominated by 
non- gravitational forces; the couplings gi between the 
various scalar fields are independent of the gravitational 
sector. Moreover, the homogeneous parts (pi, I > 1, typ- 
ically grow during the resonance phase, with the energy 
coming from the inflaton condensate, ipi. 



C. Homogeneous parts of the fields 

There is a subtlety, however, in how to define ipi for 
/ > 1 Q. One approach is to define them implicitly as 
the homogeneous objects, ipi, which satisfy the unper- 
turbed Klein-Gordon equations: 



(pi + 3H(f! +Vi=0. 



(42) 



where Vj = dV/dipj. Our main complaint with this def- 
inition is that it is not operational. It does not tell one 
how to find ipi from the full fields </>/(i, x). More impor- 
tant for the study of preheating, these equations (which 
simply express the conservation of the homogeneous part 
of the energy density) miss the gauge-invariant energy 
contributions from the k = mode of <&k- While in 
standard analyses this is unimportant, in preheating, a 
large amount of energy is transferred to this mode, which 
should be included as a homogeneous backreaction ef- 
fect. As we stressed in different terms earlier, it is im- 
portant to appreciate that because gravity has negative 
specific heat, the energy flowing into (&k=o actually aids 
the growth of field- fluctuations Sipik- 

To include this effect, we define the homogeneous parts 
of the fields J > 1 as: 



fi = S<pi,k=o- 



(43) 



Thus the background parts of the coupled fields are the 
homogeneous k = portions of the fluctuation spectra, 
bmce, as we saw m Sec. H it is the small-A; modes which 
grow the quickest for models with broad resonance, the 
f i might grow at least as quickly as other modes Sipik for 
k =/= but small. Then the combination (piSipik, which 
acts as a source for grows even more quickly than 
either term alone. This, in turn, makes fik,eS for the 
modes $^ greater than that for for any of the matter 
fields. We confirm this behavior below, in Sec. VI. In the 



however. In Sec. VIC we employ the definition implicit 
in Eq. (^|) and show that the resonances remain, if not 
quite as strong. 

Now with the definition (f43|), even if all k > modes 
of each of the fluctuations Sipik happen to lie within sta- 
bility bands at a time t, $fc will still grow if one of the 
homogeneous fields tpi lies within a resonance band. It 
is a simple combinatorial problem to see that as the to- 
tal number of fields, N, increases, the probability that 
all fields ipi and Sipjk lie simultaneously within stabil- 
ity bands, and hence that $^ receives no growing source 
terms, decreases quickly. 

As an alternative operational definition, one might use 
the field variances as a basis for defining the homogeneous 
fields (pi (I > lg viz.: 



"I 1/2 



(27T)- 3 / d'k^JkY 



(44) 



While the growth rate of (pi with this definition may be 
less than that when using Eq. (43), it is sure to exhibit 
resonance since the fc-space integral must cross all the in- 
stability bands. This is reflected in all numerical studies 
of preheating which show the rapid growth of the vari- 
ances ||. Hence the resonant growth will still exist for 
the metric perturbations 

Even using spatial averaging of the fields 4>i(x,t) to 
define the background fields ipi would lead to resonant 
growth of the tpi, since the averaging again covers insta- 
bility bands when viewed in Fourier space. ^ 

As we will see in Sec. VI, these choices for how to 



define the homogenous background fields affect quanti- 
tative details, but not the qualitative behavior, of the 
resonant amplification of modes during preheating. 

In contrast, space- time averaging, as used by Hamazaki 
and Kodama |^l[], removes the oscillations (and hence 
the singularities of their particular evolution equations) 
of the fields. This is effectively equivalent to taking the 
limit q — > and no resonant production of entropy per- 
turbations was found. Since the perfect fluid equiva- 
lence already exists with oscillating scalar fields, and non- 
singular systems exist with which to study the problem, 
we disagree with their procedure and note that this is per- 
haps the single most important reason why they did not 



**See e.g. Ql f° r a recent implementation of this definition 
in inflation. 

^Actually the situation is rather more subtle: no known 
covariant and gauge-invariant averaging scheme exists, even 
within classical general relativity. However, the above state- 
ments are true for volume averaging in synchronous coor- 
dinates and will remain true for any reasonable averaging 
scheme. See ]46| for a discussion of backreaction and aver- 
aging issues relevant to inflation. 
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find the super-Hubble resonances. Taruya and Nambu 
p2fl . on the other hand, include the inflaton oscillations 
and find resonant growth of the entropy perturbations 
(isocurvature mode). They claim however that it has no 
effect on the adiabatic mode, which they claim is con- 
stant. Further, they do not allow the resonance to go 
nonlinear and hence claim that the isocurvature mode 
decays to zero, essentially leaving no observable traces. 
Contrary to this, our explicit simulations show that 4> is 
strongly amplified and once it becomes nonlinear, which 
on general grounds occurs at or before the time at which 
the matter-field fluctuations do |0|, the role of the ex- 
pansion in damping the amplitude of the modes is not at 
all clear. 



V. CAUSALITY 

Keeping the important distinction between H' 1 and 
the particle horizon in mind p^ , we may now consider 
the role of causality in the immediate post-inflationary 
universe. Causality involves correlations across spatial 
distances. As emphasized in ]48[ | , this is best studied via 
unequal-time correlation functions in real space, rather 
than in Fourier space. In particular, causality requires 



(ST 00 (r, 77)<5T 00 (0, fj)) = for r > 77 + f, , 



(45) 



where 77 is conformal time (adrj — dt). For simplicity, 
we assume that one of the fields 4>j,J > 1, dominates 
preheating (fij > /jj, VI). The correlator ( [45| ) is then 
essentially controlled by the single field <pj. From Eqs. 
( p6| ) and (p8|), Eq. ( |45| ) shows that it is sufficient to 
study the unequal-time correlation functions of the field 
fluctuations, since each term within 5Tqq will be propor- 
tional to the field fluctuations to a good approximation, 
for modes within a resonance band. Thus, the function 
we need to study is A(r, 77,77) = {Stp(r, rj)Sip(0, fj)), i.e., 



A(r, 77,77) 



d 3 k 
(2tt) 3 



(46) 



Inside a resonance band, these modes will grow as 
^fkiv) — Pk{il)e^ kn , where Pkiv) is a quasi-periodic, 
decaying function of time. Evaluating Eq. ( [l6|) in the 
saddle-point approximation, and using Eq. (11.4.29) of 
11, yields 



A(r,7?,^) w P^(T])P k (fj)e 



exp( 



27r 3/2£3 



(47) 



where fc max is the wavenumber at which fj,/. is maximum, 
and 



£ 2 =4(77 + 77) 



d 2 Vk 



dk 2 



(48) 



Causality thus places a constraint on the characteristic 
exponent: 



dk 2 



1 

<4" 



(49) 



Note that causality restricts the shape of the spectrum 
of amplified modes, but not directly the wavelength of 
the perturbations. The approximation sometimes used 
within studies of preheating, that the distribution of 
amplified modes falls as a spike, S(k — fc r0 sonancc), vio- 
lates causality, since this requires that the field fluctu- 
ation contain correlations on all length scales, even for 
k/a ^> H. Thus, subject to Eq. ((49|), the amplification 
of long-wavelength perturbations may proceed strictly 
causally. This raises the question of whether or not met- 
ric perturbations, amplified at preheating, could affect 
observable scales even today. 

The result in Eq. ( ^) is of course limited by the single- 
field and saddle-point approximations, and also by the 
fact that the true growth in resonance bands will not 
be so cleanly exponential for multi-field cases in an ex- 
panding universe, since, as discussed in Sec. [n], realis- 
tic models of preheating have time-dependent fJ,k,es(t)- 
Thus the remaining exponential tail in Eq. ( p9| ) is an ar- 
tifact of the various approximations used. The key point 
is that, regardless of dynamical details and approxima- 
tion methods, causality is not a question which can easily 
be answered in Fourier space, and no naive estimates of 
constraints on smallest k modes which could be amplified 
causally can be made, e.g. those based on comparisons 
of k with H . 

The point we wish to emphasize is that the field equa- 
tions are relativistic, hence causality is built in."** The 



governing equations of motion, by themselves, make no 
distinctions between k/a < H and k/a > H regimes. 
A solution which satisfies the initial and boundary con- 
ditions must, therefore, be causal [see Fig. (||)]. No 
matter or energy is moved superluminally to produce 
the super-Hubble resonances at preheating. Any appar- 
ently acausal behaviour must instead be due to the ini- 
tial conditions, which, as we have emphasized above, arc 
crucial for stimulating the super-Hubble resonances we 
have discussed. The resonant amplification of very long- 
wavelength perturbations at preheating stems from the 
initial coherence of the inflaton condensate at the end 
of inflation. We return to the question of the inflaton 



condensate's coherence below, in Sec. VII 



As a final demonstration of this point, let us consider 
an inflaton distribution which is not a delta-function at 
k = but is rather distributed over fc-space. For clarity, 
let us ignore the metric perturbations completely. In this 



*We thank John Barrow for making this point to us. 
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case, the evolution equation for \k with a g<t> 2 x 2 coupling 
is of the form: 



Xk + 'iHxk + ^Xk 



9 



(27r) e 



d 3 k'd 3 k"c 



Ok-k'Xk-k' 



= 0. 



(50) 



We see immediately that if fa = for k < k clit , then 
there is no resonant growth of Xk modes for k < fc cr it 
since the inflaton is not correlated on those scales. From 
this it is easy to see that Xk modes are amplified by res- 
onance only up to the maximum scale on which the in- 
flaton is correlated. Of course, inverse cascades (see Sec. 



VII) can in principle amplify longer- wavelength modes, 
but this is another, essentially nonlinear, process. 




t = const. 



Coherently Oscillating 
Inflaton Condensate 

FIG. 3. Schematic plot illustrating how super-Hubble 
modes can be amplified without violating causality. To solve 
the horizon problem the inflaton must be correlated on vastly 
super-Hubble scales at the end of inflation, and hence this zero 
mode oscillates with a spacetime-independent phase. During 
preheating, fields at each spacetime event react only to physics 
in their causal past, which since the end of inflation may be 
very small. However, because the inflaton condensate is ex- 
actly the same over vastly larger regions than H -1 at the end 
of inflation, super-Hubble modes can patch together the infla- 
ton behavior in each region to reconstruct ("feel") the glob- 
ally oscillating condensate and hence grow resonantly. Due 
to the isometries of the spacelike surfaces of constant 0, all 
such patches are equivalent and can be mapped onto one an- 
other by a translation. Only once second-order backreaction 
effects destroy the coherence of the inflaton oscillations are 
these super-Hubble resonances weakened and perhaps ended. 

The linear case which we solve numerically is qual- 
itatively very different from the fully nonlinear regime 
where solutions are not known, as occurs for example 
in the studies of topological defects [Q. In that case, 
causality can be imposed as a constraint to determine 
whether a proposed ansatz is physical (i.e. on-shell) or 
not. This has been used to great effect in recent years in 
studies of the evolution of the energy-momentum tensor 
for causal defects, and in that context, where perturba- 



tions must come from displacement of energy, causality 
sets up strong constraints on the possible metric pertur- 
bation power spectrum. This is described in detail in 
MM- 



VI. NUMERICAL RESULTS FOR THE 2-FIELD 
SYSTEM 

We now study a two-field model involving a massive 
inflaton, (j>, and a distinct, massless scalar field, Xi with 
the potential 



(51) 



Taking x to be massless is realistic in the sense that there 
are a large number of effectively massless degrees of free- 
dom > 100) at these energies in all post-Standard 
Model theories. If g > 0, the minimum of the effective 
potential in Eq. (|5l]) is \ — 0, which we assume holds 
during inflation. For times after the end of slow-roll in 
this model, a good approximation for the evolution of the 
inflaton background field is: 



tp(t) w (p (t 



sin (mi) 



mt 



(52) 



This approximation improves with time and is quite ac- 
curate for ip < M p \. We show it in the inset to Fig. (Q). 
This ansatz for the oscillating background field appears 
to provide two free parameters: fo(to) and to. Actually, 
there exists only one real degree of freedom here, deter- 
mined by the physical amplitude of the inflaton at the 
start of preheating. We have fixed <£o(*o) = 0.3M p i, so 
that taking t Q — > (which gives sin(mi ) / mt — > 1) cor- 
responds to a standard chaotic inflation scenario [gj; we 
may think of this as a large-amplitude limit. As becomes 
clear in the following simulations, the resonances are 
more effective and modes $k go nonlinear more quickly 
for smaller to; in particular, see Fig. (|l^) for a simulation 
with mto = 1. In this case, relevant for simple chaotic 
inflation, modes go nonlinear after only mAt ~ 50 
even for quite weak values of the coupling g (and hence 
of the resonance parameter q ~ O(100)), much earlier 
than field modes 5<pi would go nonlinear in this simple 
model when $ is neglected (compare with ||). 

Increasing the arbitrary parameter to lowers the initial 
amplitude of inflaton oscillations, and hence allows us to 
study other, low energy-scale inflationary models, some 
of which we discuss further in Sec. IX. Parametrizing 



our equations this way means that the quoted values for 
the resonance parameter 



9^1 



(53) 



correspond to the values for to = 0; larger start-times, 
to > 0, will correspond to smaller q c g < q. Note that here 
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q differs by a factor 4 from the q that usually appears in 
the Mathieu equation. 

The Friedmann equation gives an approximate scale 
for the Hubble radius at the start of preheating, H(to) ~ 
2am, where m is the inflaton mass, and a gives the 
physical amplitude of the inflaton at the start of pre- 
heating: <p(to) = aM v \. For chaotic inflation, this 
yields H(t^)/m ~ 0.6 at the start of preheating; modes 
with k/a(to) < H(to) are thus, roughly speaking, super- 
Hubble modes. For our simulations, we then use the 
Hubble expansion rate as averaged over a period of the 
inflaton's oscillations, H — 2/(3t), which is appropriate 
for a massive inflaton. tn We also take the homogenous 
portion of the x ne ld to be x = (<^x)fc=o; as discussed 
above, using y/ (Sx 2 ) or the spatial average of x would 
make quantitative changes, but would not alter the gen- 
eral features or qualitative behavior presented here. 

In each of the simulations presented below, the initial 
conditions for the metric perturbations are taken to be 
the generic, scale-invariant spectrum expected to be pro- 
duced during inflation. That is, we se t $fc(to) = 10~ 5 for 
all k. As discussed below in Sec. VH, changing these ini- 
tial conditions does not, in general, alter the qualitative 
behavior presented here. In particular, attempts to fine- 
tune these initial conditions for so as to avoid nonlin- 
ear evolution during preheating are destined to fail. We 
have also used the same initial conditions for the matter- 
field fluctuations, 6ipik(to) = 10~ 5 and (S<fiik(ta))' = 0, 
where these are measured in units for which G = 1, or 
k 2 = 8ir. The initial conditions were chosen to be scale- 
invariant (i.e. independent of k) to highlight the strongly 
fc-dependent features of the resonance structure. We re- 
turn to the discussion of the physical impact of differ- 
ent choices of initial conditions on the resonances below, 
in Sec. 



VIII D 



Note that the figures in the following 
subsection plot natural logarithms (base e), rather than 
base-10. 



A. Early Growth of Perturbations 

In our numerical simulations, we integrate the coupled, 
linearized equations of motion for 8cj>k, X^ and $Xk, 
based on Eqs. ( [i"o|) and (41), with the specific poten- 
tial of Eq. (|5l|). The resulting 8-dimensional system 



The use of such time-averaging for the expansion rate and 
similar background quantities is acceptable in models domi- 
nated by non-gravitational terms, as is the case for the multi- 
field model we study here. This is in general not a good ap- 
proximation when studying gravitational waves or fields with 
strong non-minimal couplings to the Ricci curvature scalar, 
since in these cases oscillations in the Ricci curvature can lead 
to strong resonances Ml. 



of first-order ordinary differential equations were inte- 
grated using a 4th-order Runge-Kutta routine with vari- 
able step-length automatically implemented to ensure ac- 
curacy. Two equations are required each for the evo- 
lution of 6 (fik, &Xk, and Xi an d two for one with 
k > and one with k = 0. We use Planck units for 
which G = 1, k 2 = 8ir, and a dimensionless time- variable 
z = mi. Given the background expansion H — 2/(34), 
the scale factor grows on average as a oc z 2 / 3 . In other 
numerical studies of preheating, in the absence of metric 
perturbations, rescaled conformal time r = mr\ is often 



used; in this case, 



r , which facilitates comparisons 



between our results and earlier results. 
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FIG. 4. <E> vs mt for q = 5 x 10 3 , for the k = mode. Note 
the initial resonance which invalidates the linearized equations 
of motion within a few oscillations. The linearized solution 
then damps away (after mt ~ 300) due to the expansion, 
but this cannot be taken as indicative of the evolution of the 
nonlinear solution. Interestingly there are two characteristic 
frequencies of oscillation of after mt ~ 250 - one high 
frequency and one low frequency. The inset shows the evolu- 
tion of the inflaton condensate and the damping due to the 
expansion of the universe. 

Our numerical routine was tested against several lim- 
iting cases: (1) The standard results for preheating when 
metric perturbations are neglected, in which the Klein- 
Gordon equations for the fluctuations decouple at linear 
order and the system is 4-dimensional with $ set to zero. 
In this limit, the standard results of |5|||^9| are repro- 
duced, as discussed further below. (2) Limiting regions 
of the parameter space with q = were also checked. We 
neglect backreaction on the inflaton and nonlinear mode- 
mode rescattering in these simulations, but include the 
linear-order coupling which results between 5<fik and 
Sxk- The truly nonlinear effects are of course crucial for 
ending the first, explosive preheating ph ase, and we dis- 
cuss some of these effects below, in Sec. VII. 
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FIG. 5. The evolution of ln^Xfc for Q = 8000 and k = 
(top), k = 30 (middle), and k = 100 (bottom). Note the 
significant resonant growth of the k — mode. 

First consider the growth of metric perturbations in 
our model. Fig. (|l|) shows the evolution of two different 
modes, k — and k = 20, as measured in units of m. The 
first is clearly super-Hubble scaled, whereas the second 
lies well within the Hubble radius. As is clear in this 
figure, both modes evolve quasi-exponcntially for certain 
periods of time; that is, the modes slide in and out of 
various broad resonance bands. The fc = mode in fact 
begins to grow considerably after a time mAt « 40, after 
fewer than seven inflaton oscillations. 
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FIG. 6. 8<j>k vs mt evolution, q — 8000. Top: note the 
complete lack of particle production in the absence of metric 
perturbations. Lower: the particle production due to linear 
rescattering via <&k terms is clear. (Note also that for the 
inflaton field, we maintain a distinction between tp, the oscil- 
lating, coherent condensate of zero momentum, and 5</>k=o- It 
is the latter which we plot here.) 



Figures (|l|) and (|12|) reveal the strongly scale- 
dependent spectrum for which results from the early, 
linear-regime evolution of the coupled system. This is the 
analog for $^ of the strongly non-thermal distribution of 
particles produced in the first stages of preheating. We 
will return to the question of the power spectrum for $fe 
during the preheating epoch in Sec. 



VII 



Of course, the 

evolution for both modes plotted here beyond (t) ~ l 
cannot be trusted based on our linearized equations of 
motion, and nonlinearities will end the resonant growth. 
We plot the strong growth here as an indication that 
Floquet theory remains an effective tool for describing 
the linear regime. Rather than the simple Minkowski 
spacetime results, the effective Floquet indices Hk,eft(t) 
become time-dependent: the modes $^ slide into and 
out of resonance bands, each characterized by a certain 
exponential-amplification rate. 

The coupled scalar field fluctuations 8\k also experi- 
ence dramatic quasi-exponential amplification when the 
coupling to $ is taken into account, as shown in Figs. 
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We may now compare the behavior of the matter-field 
fluctuations, 5<f>k and 5\k, when we include the coupled 
metric perturbations, and when we neglect these pertur- 
bations as in the standard preheating literature. Fig. 
(Q) contains a plot of the inflaton fluctuations 5<pk with 
and without including the coupling to <E>/c, and Fig. (Q) 
contains a similar plot for 8\k- 
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FIG. 7. Here we compare the evolution of 8\k vs mt (k — 0, 
q = 8000, mto = 30) with and without the full metric cou- 
pling. 

Note that for the early times plotted here, the behavior 
of each of these fluctuations when we neglect matches 
well the behavior found in |]J| . To facilitate comparison, 
recall that z oc r 3 , and that we plot here \ — (8x)k=o 
rather than (Sx 2 )- From Fig. (|5|), it is clear that the 
k = mode of 8xk dominates the growth, and hence 
contributes most to (Sx 2 ), so this makes for a fair com- 
parison. The gradual decay of both 5<fik and Sxk shown 
here for early times, when <E> is neglected, matches the be- 
havior found in |5|||. In these earlier studies, significant 
amplification of each of these fluctuation fields occurred 
only for r > 10, or, roughly, z > 10 3 . Instead, when the 
coupling to the metric perturbations is included, we find 
significant growth in each of the field fluctuations well 
before these late times. 

The strong, resonant growth of both S<fik and 5xk con- 



firms that the metric perturbations serve as both a source 
and a pump for the matter-field fluctuations. We con- 
jectured this above and in Paper I, based both on the 
form of the perturbed energy density and on the linear- 
order gravitational rescattering. This is a significant new 
effect which simply cannot be ignored when studying re- 
alistic models of preheating, which involve more than one 
scalar field. Broad-resonance amplification persists in an 
expanding universe when q c s ^> 1 , and this amplification 
now affects the fast-growing <f>fc modes then further 
stimulate growth in 5<f>k and b~Xk- 

Moreover, note that the rate of growth Hk,eS for $fc = o 
is greater than that for both 5<fik and for 8xk- This, too, 
is consistent with our analysis above, based on the form 
of the constraint equation, Eq. (^l|), for 

Again, we emphasize that the point of these figures is 
not to study quantitative details of the field and metric 
perturbations at late times, which is impossible using our 
linearized evolution equations. Rather it is that quasi- 
exponential growth of both of these types of perturbations 
will take place early after the start of preheating, and 
long before it takes place when the metric perturbations 
are neglected. 



B. Time to nonlinearity 

Given the rapid, quasi-exponential growth of these per- 
turbations, the most appropriate physical quantity to 
study is the time t n \ it takes a mode k to saturate its 
linear theory bound, that is, |3>fc(i n i)| = 1- Clearly in a 
system dominated by gravitation alone, t n i ~ Hq 1 , the 
current Hubble time, which is much longer than reheating 
can last. At preheating, however, in multi-field models 
in the broad-resonance regime, we expect instead that 



t 



ni 



(54) 



Modes lying within a stable band have [ik = and t n \ — 

00, whereas modes within a resonance band will saturate 
their linear-theory limit at some finite time. 

In Fig. (H) we show t n \ versus the resonance parame- 
ter q for moderate values of q and with a cut-off in the 
integration time of to At = 250. There is a clear phase 
transition at q ~ 2800, where modes start going nonlinear 
for the first time. This is followed by a succession of nar- 
row stability bands and a second phase transition around 
g* » 6000 where no more stable modes are found. For 
q > q*, essentially all modes go nonlinear by to At ~ 100, 

1. e. after only a few inflaton oscillations. 
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FIG. 8. Time to nonlinearity t n i vs q for the k — mode of 
<3?fc . Note the transitions at q ~ 2800 and g, ~ 6000. The inset 
shows t n l vs g for negative q and fc = 0. There are no stable 
bands and the time to nonlinearity is short. The numerical 
solution for the negative-q case is closely approximated by the 
simple scaling law t n \ ~ 1/ \/\q\ predicted by Eq. (|5f]) . 

The inset figure shows t n \ for the case of a negative- 
coupling instability amongst the field modes, which pro- 
duces a more efficient resonant amplification of all of 
the perturbations. Note in particular that t n \ matches 
very well the prediction from earlier studies of negative- 
coupling instabilities that (29[ 



-1/2 



(55) 



In Figs. (|9f]lC|) we show the effect of decreasing the 
starting value of to. Since decreasing t increases the ef- 
fective inflaton oscillation amplitude, we expect that it 
should increase the power of the resonance and hence de- 
crease t n \ and the number of modes lying in stable bands. 
This is clearly seen in the figures. Note in particular 
that for most values of q in the range 1500 < q < 5000, 
and with 5 < mto < 30, modes $fc go nonlinear after 
mAt ~ 40 — 100, well before any such nonlinearities ap- 
pear in studies of preheating which neglect the coupling 
to $. 

Note further that as t is increased, decreasing the 
strength of the resonance, the resonance band structure 
becomes more and more clear. This behavior of t n \ re- 
veals important qualitative agreement with earlier nu- 
merical studies of preheating. For preheating into mass- 
less fields x, significant amplification in an expanding 
universe had only been found for q > 10 3 j5). The rapid 
growth of $k modes found here for our model of the cou- 
pled, massless 'fields' $ and Xi f aus i n this range for q. 
Moreover, as q e s increases, the dependence of t n \ on q 
becomes more fully stochastic, as seen especially in Fig. 
(pd|). In this figure, no stable bands persist for q > 3000, 
and t n \ varies non-monotonically with q. These features 
again match the qualitative behavior found in (^] for the 
stochastic, broad-resonance regime in multi-field models 



of preheating that neglect metric perturbations. We find 
this same kind of behavior here for the evolution of 
with q 1. 
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FIG. 9. tnl vs q for the k mode and for to = 10 (main 
figure) and to = 20 (inset top right). From this it is clear that 
varying to is not exactly the same as simply varying q essen- 
tially due to the 4> terms which are large when to is small (the 
region where the amplitude of oscillations drops very rapidly) . 
The large q, large to regime is different from the smaller q, 
smaller to regime which is closer in spirit to the stochastic res- 
onance regime. This is evident between this figure and figure 
([lo|) which shows ordered resonance and stable bands with 
mto = 30. 
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FIG. 10. tni vs q for k = and mto = 30. Note the 
ordered resonance bands. Inset: t n i vs q for mto = 5. At 
this small value of mto modes go nonlinear rapidly and no 
resonance band structure is apparent. 

A similar relationship governs the dependence of t n \ 
on k, as shown in Fig. (|l^). Again, note the absence 
of stable bands in /c-space for this q ^S> 1 regime. Note, 
too, that the most rapid amplification, and hence the 
smallest t a i, occur in the k — > limit. Moreover, this 
rapid amplification persists smoothly for modes greater 
than and less than the Hubble scale, with t n \ varying 
little in the range < k < 15. This provides a further 
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indication that H 1 is not a fundamental quantity for 
preheating with more than one matter field. 
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FIG. 11. t n i over a wide range of q > with mto = 30. 
The essentially stochastic variation of t n i is clear, as is the 
transition to nonlinear behaviour for q ~ 3000. t n i for q < 
is not included here. 
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FIG. 12. tni vs fc for q = 8 x 10 3 (top left inset), g = 10 4 
(main figure) and q — 7 x 10 4 (bottom right inset). It is clear 
that the small k modes are, on average, the most strongly 
amplified (with no violation of causality). This is due to 
their low momentum and hence ease of violating the adia- 
batic condition, Eq. (|ll]). Note that, as expected on general 
grounds from Floquet theory, there exist (q, fc)-parameter re- 
gions where increasing q leads to an increase in t n i showing 
that the resonance is not monotonically increasing in strength 
with q. 



C. Numerical results with an alternative definition 
of the background x 

A reasonable concern the reader might have is whether 
the strong super-Hubble resonances we found in the ear- 
lier section are an artifact of our definition of the homo- 
geneous component of x as the k = mode of 8\k- As 
already explained, we made that choice so as to be sensi- 



tive to the considerable energy in $>k=o, and as such may 
be regarded informally as an attempt to include some 
backreaction effects. 

The standard approach is to require that the back- 
ground fields satisfy the unperturbed Klein-Gordon equa- 
tion, which simply expresses energy conservation at zero 
order. In that case, and for the potential in Eq. (p3l|), we 
have: 



X + 3H X 



g4> 2 x 



o. 



(56) 



which misses the driving terms involving the homoge- 
neous $fc = o mode. The k = mode is not in itself 
gauge-invariant j5jj, but is physical, since, in the longi- 
tudinal gauge which we use, all gauge freedoms are fixed 
Q . From a physical point of view, the background quan- 
tities would correspond to the coherent condensate, with 
k ~ 10~ 28 (see discussion at the end of Sec. QI), and 
hence are strictly gauge-invariant. It is clearly important 
to demonstrate that the super-Hubble resonances persist 
for the above, standard, definition of \- I n Fig- ( B ) 
we show the evolution of $fc and 5(f>k for k — 1CP 6 for 
q = 8000, t = 10. In Fig. @ we show the evolution of 
5\k and the homogeneous X found from Eq. ( |56| ) for the 
same parameter values. Note that although X grows rel- 
atively slowly, the resonances in the k = 10~ 6 modes are 
still strong, with each of the perturbed quantities going 
nonlinear after mAt ~ 65. 

It should be clear that none of the essential predictions 
have changed, as we anticipated earlier in our discussion 
of the constraint equation. The g(j) 2 X term in Eq. ( |56"| ) at 
large q is sufficient to ensure resonant growth of X , which 
leads to resonant growth of . 

In Fig. ( |l5| ) we show t a \ vs q for k = for mto = h 
where X satisfies Eq. ([56]). Note the fine-structure and 
the sudden transition to nonlinearity, characteristic of 
the earlier t n \ vs q figures. In the usual chaotic inflation 
scenario, with mto = 1, modes <&k go nonlinear after 
mAt < 50 even for low-g values of q > 200, that is, even 
for very weak couplings g > 10 ~ 8 . 
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FIG. 13. S Xk vs rat for q = 8000 and k = 10~ 6 and 
mto = 10. Inset: The evoluion of the homogeneous mode x 
vs mt for the same parameter values and using the evolution 
eq. (p6|). Notice that although \ remains small throughout 
the simulation, all the perturbations Sxk, 5cj>k and have 
still gone nonlinear. This demonstrates that the super-Hubble 
resonances are robust predictions and not sensitive to the def- 
inition of the homogeneous fields, though the precise times to 
nonlinearity are. 

In conclusion, the existence of resonances for $fc is not 
sensitive to the choice of homogeneous mode for \- Sat- 
uration of linear theory and the need for nonlinear study 
appears to be a robust prediction of strong preheating, 
independent of the definition of the background fields, as 
long as they are sensitive to the non-gravitational forces 
implicit in preheating, as indeed they must be. 




mt 

FIG. 14. $ fc vs mt for q = 8000 and k = 10" B and 
mto = 10. Inset: The evoluion of 8<f>k vs mt for the same 
parameter values. 
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FIG. 15. t n \ vs q for mto — 1 and k — using the conser- 
vative definition of \ which satisfies Eq. (^). The complex 
band structure is still visible. 



VII. BACKREACTION, NONLINEAR ASPECTS 
AND IMPLICATIONS 



A. Nonlinearity vs fine-tuning 

In the standard theory based on conserved quantities 
such as £, the values of <E> during and after inflation are 
constrained by a simple relation such as |53] : 



1 



1 



Wf 



(57) 



where w = p/p. Clearly, 3>(t/) ^> $(U) if Wi ~ — 1, as 
occurs in slow-roll inflation. By fine-tuning <p towards 
zero during inflation, the amplification of $ can be made 
arbitrarily large. In the standard theory, this doesn't 
imply that perturbations have gone nonlinear. Rather, 
the problem is run backwards: the $ ~ 10~ 5 required by 
the CMB is used to determine the value of $(tj) during 
inflation, which therefore puts constraints on the value 
of 4> during inflation. No matter what the amount of 
amplification, the perturbations can be kept linear by 
sufficient fine-tuning of the intitial value, 

Have the previous sections simply been a discovery 
of a more elegant, realistic, and refined version of this 
amplification? And why are we claiming that the same 
fine-tuning cannot be applied in the realistic case of pre- 
heating? That is, despite the extra amplification due to 
the entropy perturbations (which invalidates the use of 
conserved quantities such as £), why can't all the reso- 
nant amplification be swept under the carpet, and $ kept 
linear, by extra fine-tuning? 

Clearly this is a crucial point. The extra amplification 
due to the quasi-exponential growth of entropy perturba- 
tions means that additional fine-tuning would be needed. 
This would exacerbate the already unpleasant fine-tuning 
of the couplings in many inflationary models, particularly 
simple chaotic inflation models. 

Yet fine-tuning, while not desirable, is however an is- 
sue of taste. Much more serious in any attempt to stop 
the perturbations from going nonlinear are the following 
points: 

(i) First, fine-tuning the value of &(ti) downwards dur- 
ing inflation will typically require making the couplings 
(to 2 , A) of the infiaton field smaller. The dependence 
of &(ti) on these parameters is typically power-law, e.g. 
AT/T tx m/Mpi in the quadratic potential model, and 
AT/T on A 1 / 2 in the quartic potential model. However, 
assuming (as is true in simple chaotic inflationary mod- 
els) that these parameters also control the oscillation fre- 
quency of the infiaton during reheating, we see that the 
resonance parameter q cx m~g is significantly increased 
when the effective mass is decreased. This makes the res- 
onance stronger and increases the effective Floquet index 
Mfc,cff ■ Since depends exponentially on [ik,cS, it is un- 
likely that any fine-tuning gains made by reducing the 
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coupling constants of <f> will survive through the end of 
the resonance in reheating. 

(ii) Second, and perhaps of more fundamental impor- 
tance, concerns the timing of the end (tend) of the ini- 
tial resonance phase. Floquet theory guarantees that the 
fluctuations will grow resonantly until the end of this 
linear phase. Typically in preheating studies tend 
has been taken to be the time when backreaction be- 
came 'important,' which was usually taken on dimen- 
sional grounds to be when </> 2 ~ (x 2 ), (S(f> 2 ), whichever 
came first. In Paper I we showed that when the metric 
perturbations are included in the analysis, one expects 
backreaction to become important in the field fluctua- 
tions and metric perturbations around the same time. 
Our simulations here suggest that at large q (in the multi- 
field case), it is often the metric fluctuations <I> which go 
nonlinear first. 

In the absence of special circumstances (such as a 
strong x self- interaction, A x ), the end of preheating will 
occur when second- and higher-order nonlinear processes 
dominate the linear resonance, no matter what degree of 
initial fine-tuning is employed. Thus, regardless of spe- 
cial fine-tuning during inflation, preheating in multi-field 
models in general will proceed at least until <£> enters the 
nonlinear regime. For this reason, we consider some im- 
plications of this nonlinearity in the remainder of this 
section. 

B. Mode- mode coupling and turbulence 

Several backreaction effects have been highlighted in 
the preheating literature. Two of them may be approx- 
imated with Hartree-Fock or large-iV schemes: (i) the 
change to the inflaton amplitude of oscillation, ipo, and 
(ii) the change to the inflaton effective mass, m c s = m 2 + 
g 2 (x 2 ), and hence to the frequency of inflaton oscillations 
lUpPSfl . The first is a strongly non-equilibrium effect, be- 
ing governed by non-Markovian dynamics rather than by 
the simple Markovian approximation of adding a T(f> term 
to the inflaton equation of motion M. In both of these 
cases, the effects derive from the growth of the backreac- 
tion term S introduced in Sec. O, S x ~ J d 3 k\xk\ 2 and 
^<t> ~ / d 3 k\5(f)k\ 2 - The implications of both (i) and (ii) 
for the resonance parameter q appear straightforward: 
given q oc ip^/m 2 , the resonance parameter will fall as 
the inflaton amplitude falls; it will likewise fall as the in- 
flaton's effective mass grows. If the effective mass of the 
inflaton should decrease, then these two effects would 
compete. 

As treated within the Hartree-Fock or large- N approx- 
imations, these backreaction effects are essentially inde- 
pendent of fc; given the form of £/, they simply sum up 
the growth of modes within certain resonance bands, as a 
quick way to track the transfer of energy from the coher- 
ent inflaton into various k modes. As discovered in lattice 



simulations Q and explained analytically ||, these ap- 
proximation schemes thereby miss crucial elements of the 
nonlinear evolution. In particular, fc-dependent mode- 
mode coupling can play a large role before the end of 
preheating. It is clear from experience that any nonlin- 
ear field theory may induce mode-mode coupling, and 
this is the origin of rescattering in ordinary preheating. 

In a Hubert space, the Fourier transform of a product 
is simply the convolution of the individual Fourier trans- 
forms; i.e., T{f ■ g) — T{f) * F{g). The conversion of a 
nonlinear term in the field equations into Fourier space 
therefore yields one or more convolutions of the spectra 
of the fields involved. The nonlinear term g<fi 2 x m the 
equation of motion for for example, becomes the dou- 
ble convolution in Fourier space (which is symmetric in 
the arguments of <fi and x): 

g J J d 3 k'd 3 k"</> k _ k /_ k »0 k _ k /Xk- 

Assuming statistical isotropy and homogeneity, this re- 
duces to a double scalar integral with measure k' 2 dk' ■ 
k" 2 dk". It is precisely this coupling of modes of differ- 
ent momenta which is missed by the Hartree-Fock and 
large- TV approximation schemes. 

In the case of perturbations around FRW, we must 
make a splitting of the scalar matter fields into homoge- 
neous background fields and fluctuations. When this is 
done, two kinds of terms arise in performing these con- 
volutions: terms for which both k — k' — k" and k — k! are 
non-zero (i.e. scatterings involving two non-condensate 
particles) , and terms in which one of the <f>k particles be- 
longs to the condensate, with k = 0. The latter type of 
term is dominant at the early stages of mode-mode cou- 
pling, because the inflaton condensate occupation num- 
ber is still high H,|. 

When metric perturbations are included, similar con- 
volution terms (and hence rescatterings) appear beyond 
linear order, now involving the modes of in addition 
to Sip ik- At first glance, it might appear simply too diffi- 
cult to study these terms analytically, since at second or- 
der, the equations describing the coupled system become 
tremendously complicated. The first complication arises 
because the variables we use, such as $ and Sip (and in- 
deed any other known variables) , are no longer gauge in- 
variant. Thus a gauge-dependent formalism must be used 
in practice, although in principle some gauge-invariant 
generalizations of these variables might be found p3[ . 
Further complications arise from the non-zero pressure 
and the increased variety of cross-terms among the many 
scalar matter fields. These issues in gravitational per- 
turbation theory have been explored in a cosmological 
setting in studies of dust Einstein-de Sitter spacetime 
II- 

These complications, of course, are bound to stymie 
any analytic progress on the fully nonlinear evolution of 
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multi-field systems at preheating. Yet despite these com- 
plications, we may make a survey of the kinds of terms 
which will be relevant to the specific question of mode- 
mode coupling and rescattering at second order. At sec- 
ond order, the system will contain convolutions of two 
first-order terms, in addition to new terms formed from 
background quantities multiplied by explicitly second- 
order objects. Because the background quantities yield 
delta-functions in Fourier space, however, these latter 
terms make no contributions to mode-mode coupling. 
Thus, modulo subtleties regarding gauge fixing, the non- 
linear effects from mode-mode coupling arise from com- 
binations of terms which are now familiar from the lin- 
earized equations studied above. 

Consider, as an example, possible effects of mode-mode 
coupling for a single- field model with small-amplitude os- 
cillations Upo M p i). In this case, the linear theory 
predicts [[L8| no amplification of fc — > modes of $ at 
preheating. The constraint equation at second order for 
this simple case would take the form: 



(2) 



2T*« = 



dk' k' 



x-Wx (!) 



(58) 



The integral on the right-hand side can transfer power 
from high-fc to small- k modes (an inverse cascade), as 
well as from small- A: to larger k (a forward cascade). Such 
inverse cascades could in principle amplify super-Hubble 
modes of $ even in such a small-amplitude, single-field 
model, based on these mode-mode couplings. Inverse 
cascades of precisely this type are common in nonlinear 
systems Ipq] , and indeed evidence for them was found 
recently in a numerical study of preheating Pq| . 

While still a topic of significant debate in turbulence 
circles, the classical results of Kolmogorov on scaling in 
turbulence (56) suggest a fc~ 5 / 3 scaling law for the energy 
cascade in 3D turbulence. Similar direct cascade behav- 
ior has been found in lattice preheating simulations ||. 
An even more radical possibility is that of inverse cas- 
cades, in which power is transferred from small to large 
scales. Such behavior is believed to occur in energy trans- 
port in 2D hydrodynamics, as well as in particle number 
in the nonlinear Schrodinger equation, which signals Bose 
condensation. In addition, there is the claim [|55| that 
energy cascades are inverse in nature in turbulence and 
magneto- hydrodynamics with power-law spectra oc fc Q , 
if a > —3, as is the case in inflation, although see also 
f37f . Of course, in the cosmological context causality 
must play a crucial role. Here, the coherence of the in- 
flaton condensate is irrelevant and cascades must obey 
causality, yielding the usual result that the part of the 
spectrum added by mode-mode coupling can go on large 
scales at best like fc 4 058). 



Similar terms involving convolutions over <4 1} and 
6ip\]^ would appear at second order in the equations of 



motion for both the gravitational and matter-field fluctu- 
ations in multi-field models. For the matter-field fluctu- 
ations, these nonlinear $ * Sipi convolutions would there- 
fore add a new source for rescattering to those consid- 
ered in earlier studies: scattering of 4> and X particles 
directly off the gravitational potential. This nonlinear 
mode-mode gravitational rescattering would complement 
the gravitational rescattering we found already at linear 
order, in the form of the VijSip^ cross-terms examined 
in Sec. IV, ultimately leading to larger variances in the 
matter-field fluctuations than the linear theory alone will 
produce. 



C. The post-preheating power spectrum and the 
CMB 

Perhaps more important, we may expect the results 
found in earlier numerical studies of preheating, regard- 
ing turbulence and the final shape of the power-spectrum 
for fluctuations H , to carry over to the metric perturba- 
tions <!> as well. As found there for the case of coupled 
matter-field fluctuations, turbulence and cascades will 
lead, in general, to a re-establishment of scale-invariance 
for the power spectra (P/(fc)) for each 8<pi. In the case 
of then, these mode- mode coupling effects ultimately 
would lessen the sharply scale-dependent spectrum which 
results at the end of the linear era of preheating, though 
as discussed in the previous subsection, causality will 
limit its effects for very small k. This nonlinear effect on 
P$(fc) could thus produce a nearly scale-invariant spec- 
trum by the end of preheating above a threshold value 
of fc. It is crucial to recognize, however, that this "late- 
time" (end of preheating) power spectrum -P$(fc) would 
be virtually independent of the scale-invariant spectrum 
produced during inflation (|59f . The nearly-flat spectrum 
measured by COBE, which probes P$(fc) over very small 
fc, would therefore be a complicated combination of these 
two spectra. 

However, based on earlier work on second-order pertur- 
bations and the CMB |30| , it is clear that amplification of 
super-Hubble modes to the stage where second-order per- 
turbations cannot be neglected is not compatible in naive 
models with the COBE DMR results that AT/T ~ 1CT 5 . 
Hence some mechanism is required to damp the spectra 
or limit the resonances before they can be observationally 
compatible. Some implementations of this requirement 
are discussed in Sec. 



VIII. 



An interesting possibility also remains for observa- 
tional consequences of P$(fc) as it evolves through its 
highly nonlinear mode-mode coupling regime in the later 
stages of preheating. This second effect concerns larger- 
fc (shorter-scale) perturbations. The turbulent, nonlin- 
ear evolution of modes at preheating could effec- 
tively mimic generic features of the power spectra pre- 
dicted by topological defect models. In particular, the 
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broad-resonance, stochastic amplification to nonlinearity 
of certain modes could potentially smear out the sec- 
ondary acoustic Doppler peaks in the I > 100 portion 
of the CMB spectrum, if the modes remained nonlin- 
ear all the way down to decoupling. These sharp sec- 
ondary Doppler peaks are expected in nearly all inflation- 
ary models (when one ignores completely the behavior of 
$ during preheating), yet are usually absent from de- 
fect models (which feature "active," incoherent spectra, 
especially in these higher multipoles) ]6lj ]. 

Combined with the spectrum at small k, this could 
produce a hybrid spectrum for P$(fc): (nearly) scale- 
invariant at scales probed by COBE (and hence naively 
similar to standard inflationary predictions, though phys- 
ically of radically different origin), yet lacking any clear 
secondary Doppler peaks at shorter scales (and hence 
naively similar to standard predictions from defects mod- 
els, though again for very different underlying physical 
reasons). The specific effects of the resonant amplifica- 
tion of $fe modes at preheating on the CMB requires 
careful study, and is currently under investigation by the 
authors. 

D. Coherence of the inflaton zero mode 

As we discusssed in Sec. [y], the super-Hubble reso- 
nances which form the basic element of this paper owe 
their existence to the coherent oscillations of the infla- 
ton zero mode. These oscillations lead to stimulated ef- 
fects that are missing if one treats the inflaton as an 
incoherent fluid ||. Hence, in addition to the backre- 
action mechanisms outlined above, we expect that an- 
other important backreaction event will be the reduc- 
tion in the correlation length of the inflaton condensate. 
This is expected on very general grounds based on the 
fluctuation-dissipation theorem |p2[ : backreaction will 
cause the long-range coherence of the inflaton condensate 
to be broken. Eventually the phase of the inflaton's os- 
cillations will become spacetime dependent, with its own 
correlation length that is expected to decrease rapidly 
& 

As a result we expect the treatment of the inflaton 
as an effectively infinite, coherent condensate to break 
down on a timescale t C oh- However, the reader should 
note that at first order in perturbation theory this effect 
does not appear, and only at second and higher order 
can the loss of coherence appear consistently. Hence the 
super-Hubble resonances cannot be eliminated. These 
higher-order effects, rather, would decrease the time dur- 
ing which the very-long wavelength metric perturbations 
would be amplified, and hence could limit the final am- 
plitude of these modes. 

For times t > £ co h, as the inflaton correlation length de- 
creases rapidly and the inflaton begins to lose coherence 
as higher-momentum modes get amplified, we expect the 



resonance for the k/a <C H modes to weaken signifi- 
cantly. This would likely end further resonant growth 
of the super-Hubble modes. From that time on, the only 
mechanism which could amplify the k = mode would 
be inverse cascades. 

The coherence issue has been discussed briefly previ- 
ously. In |Q it is argued that the inflaton's coherence is 
only lost if the decay rate of the condensate, consid- 
ered as a collection of zero-momentum bosons, is greater 
than H. In this sense, the expansion of the universe 
has a 'healing' effect on the coherence. For a coupling 
of the form in Eq. (@), they argue (using only the high- 
frequency limit to calculate the cross-section for T^) that 
this ratio is 

^ « 1.2 x 10V, (59) 

so that only for g > 3 x 10 -3 is coherence lost. Note 
that this does not depend immediately on ipo and m (m 
was cancelled in the cross-section and the average H was 
used, so ipo did not appear). If this is correct, then in 
the simplest models of chaotic inflation, Eq. (^|) predicts 
that all values of q < O(10 7 ) lead to coherent inflaton 
oscillations, regardless of backreaction. In this case, the 
super-Hubble resonances of the k sa modes would con- 
tinue until q was significantly reduced via the other q- 
reducing backreaction mechanisms. This point was con- 
sidered separately in j65| , in which it was argued based 
on semi-classical quantum gravity that coherence would 
be lost rapidly after only a few oscillations, though no 
detailed arguments were given. 

Note that while the break-up of the coherence of the 
inflaton condensate may stunt the super-Hubble reso- 
nances, it will not stop resonances on smaller scales. Fur- 
ther, recent work (see the second ref. in Qj) has shown 
that parametric resonances persist even in the case where 
the effective mass is modulated by spatially inhomoge- 
neous white noise. They further conjecture that the noise 
increases the Floquet exponent, fj,k, for all modes. If this 
is correct, the break-up of the condensate may not aid in 
ending the super-Hubble resonances at all. 

Quite apart from these gravitational considerations, 
yet another fundamental difference appears relevant be- 
tween single-field and multi-field models: the break-up 
of the super-Hubble coherence of the inflaton conden- 
sate may proceed more quickly in multi-field models than 
in the single-field case. This stems not only from the 
larger couplings between fluctuations in the multi-field 
case, but also from the greater number of directions in 
field space along which the coupled system could exit 
inflation. In other words, the inflaton condensate co- 
herence length even at the start of preheating is likely 
to be shorter in multi-field models than in the single- 
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field case.^ Similar effects have been highlighted in re- 
cent studies of tunneling in multi-field models of open in- 
flation, in which the multidimensional field-space yields 
large "quasi-open" domains inside the nucleated bubble, 
rather than a perfectly homogenous background. p6[ The 
coherence of the inflaton condensate thus remains a cru- 
cial quantity in need of further study; the amplification 
of super- Hubble modes depends directly upon it |63] . 



VIII. POSSIBLE ESCAPE ROUTES FROM 
NONLINEARITY 

Given the severity of the consequences of the super- 
Hubblc resonances with regard to observational compat- 
ibility, primarily with the CMB, a natural question, and 
one to which we have already alluded previously in the 
paper, is "can one have large q but still escape the super- 
Hubble resonances?" 



this low-energy scale corresponds to a lifetime several or- 
ders of magnitude larger than the Planck time, the "nat- 
ural" time-scale for recollapsc [jl| . This is a generic prob- 
lem of 'low-energy' inflationary models. A first phase of 
chaotic inflation which sets up appropriate initial condi- 
tions for a second phase of inflation, which in turn pro- 
duces the required density fluctuations, is an alternative 
& 

Of course this is not the only solution. The 
above objection cannot be made of compactified multi- 
dimensional models in which 4-dimensional unification 
occurs around the GUT scale, since then GUT scale in- 
flation occurs immediately upon exit from the Planck 
epoch. In the most extreme version of these models, in 
which unification occurs around the TeV scale f73|] , it is, 
however, not clear at present whether inflation can be 
implemented naturally or not p3J73]. 



B. Escape route 2: fermionic and instant preheating 



A. Escape route Is secondary phases of inflation 

The no-hair theorem [|67f assures us that should a sec- 
ond phase of inflation follow preheating, the large super- 
Hubble resonances will be smoothed with time, the uni- 
verse re-approaching the de Sitter geometry. 

Such secondary phases of inflation are desirable and 
even expected in certain models. Consider thermal in- 
flation for example |59|, proposed as a solution to the 
moduli problem with a relatively small number (~ 10) 
of e- folds. Such a phase of inflation following preheating 
would stretch and smooth, but not eliminate, the signa- 
ture of the nonlinearity achieved during preheating. The 
challenge is to identify the fingerprint that one would 
expect on the large-angle CMB in this case. 

Secondary phases of inflation are expected in multiple 
inflationary models based on supersymmetry and super- 
gravity [fzofl - Studies have already begun examining the 
effect of the sudden change in the effective mass of the 
inflaton in the transition from one phase to another on 
the CMB, though until now the change has been treated 
neglecting any resonance effects J7^j7l[] . 

Finally, it appears that the COBE results prefer a 
relatively low energy scale for the production of the 
anisotropies, U 1 /4/ £ i/4 ^ g 7 x 1Q i6 GeVj where £ = 

\M 2 X {V /V) 2 <C 1 is the slow-roll parameter (see Liddle 
and Lyth in 0). Such a low energy scale is a prob- 
lem if one also wishes inflation to truly solve the hori- 
zon/flatness problem of the standard cosmology, since 
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Until now we have considered a model of the post- 
inflationary universe based solely on scalar fields. The 
case of preheating to gauge bosons is known to be quali- 
tatively similar to the scalar case j7o| , though more com- 
plex. Fermionic preheating, on the other hand, is effec- 
tively crippled by the Pauli exclusion principle, which 
places the limit < 1/2 on the occupation number, as 
opposed to the 1 possible in the scalar case [[77| f79| . 

A natural question then is how the metric reso- 
nances are affected when the inflaton can decay only 
into fermions at the end of inflation. Our aim here is 
to present evidence that the exclusion principle effec- 
tively stiffles the resonant amplification of metric per- 
turbations at preheating, and thereby possibly offers a 
way out from the super-Hubble resonances that occur in 
the scalar case. 

Consider a massless spin i fermion, ip, coupled to the 
inflaton via a hipfiip interaction term. -0 satisfies the 
Dirac equation 



ft0(t)ty = O, 



(60) 



where 7^ are the Dirac matrices satisfying the Clifford 
algebra relations {7 /1 ,7 K } = 2n^l. Making the anzatz 
ijj = [ry^V ft +h(f>(t)]X (see @||), X satisfies the Klein- 
Gordon equation with a complex effective mass [cf. Eqs. 
(Ill) and (0)1: 



til = ul + qf - iy/qf 



(61) 



where / = (t> 2 /fo as in Sec. |nj and q oc h 2 /m 2 for a 
massive but not self-interacting inflaton, while q cx h 2 / X 
in the case of a massless, quartically self-coupled inflaton. 
The comoving number of created fermions is |79|| : 
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n>k 2 



2u>llm(X k X* k ) 



2 £7 1 



with energy density 



Pi, = / dkk 2 n k n k . 



(62) 



(63) 



Note that most of this energy density is in the form of 
inhomogeneity since nk=o < 1/2 like all other modes, un- 
less a (chiral-symmetry breaking) condensate forms with 
(ipip) 0. Assuming that such a condensate has not 
formed, one expects from Eq. ( |63| ) that the total energy 
radiated into the fermion field will be small, essentially 
due to the low occupation numbers forced by the exclu- 
sion principle. 

This is borne out by a full closed-time-path analy- 
sis of the problem |j77|, which shows that the ensemble- 
averaged energy dissipated to tp from a massive inflaton 



m h <pg 
64tt 2 



(64) 



which is rather small. This should be compared with the 
energy density in scalar \ particles at i e ndj the end of the 
linear resonance, which one can estimate in the Hartree 
approximation Q as p x (i C nd) ~ sl^l^x ~ ra 2 ip^. This 
second energy density is much larger than its fermionic 
counterpart unless h 2 > 1, and this is not allowed by 
the CMB, unless the theory is supersymmctric. Further, 
the energy pumped into the \ field predominantly goes 
into super-Hubble modes, whereas in the fermionic case, 
only a small fraction of the total fermionic energy can go 
into amplifying super-Hubble modes due to the exclusion 
principle. This is exacerbated by the small phase-space 
volume, oc fc 3 , at small k. 

The equation of state for tp must satisfy < p^p, 
and hence we are assured that 8T^ V is small relative to 
the energy density in the background. The perturbed 
Einstein field equations B 



V;($ + if$) (x6T° 



(65) 



suggest that it is unlikely that any significant amplifica- 
tion of super-Hubble metric perturbations exists in the 
case of preheating only into fermions. 

More interesting is the realistic case where scalar and 
fermionic fields coexist, coupled to the inflaton. Ex- 
amples of fermion-stimulated decay of the inflaton into 
bosons exist Jf8| , and it is not clear how metric pertur- 
bations evolve in this case. 

A particular example is the recent model known as 
instant preheating J80| , in which the sequence (j> — > \ — > 
i\) occurs, with the end point being small numbers of very 
massive fermions, with masses m x > 10 16 GeV possible. 
While it is not our aim to study the evolution of $ in 
this model, we expect that the resonances discussed here 



will not exist in their current form in instant preheating, 
since the x bosons decay almost immediately (within one 
oscillation) in the simplest scenarios, while the resulting 
fermions are few in number, and the energy transferred 
to them is typically rather small. Still, this energy may 
grow to dominate the total energy density of the universe 
with time if the inflaton energy redshifts relativistically, 
cx a~ 4 . It is difficult to say anything more quantitative, 
and these interesting issues are left for future work. 



C. Escape route 3: % mass and self-interaction 

In our numerical analysis, we limited ourselves to a 
model of preheating with a massless scalar field x- Two 
X characteristics one might expect of x are a non-zero 
mass (m x ^ 0), and self-interaction (A x ^ 0). For the 
potential 

V(<f>, x) = \™l4> 2 + h94?X 2 + \rn\x 2 + ±A xX 4 • (66) 

the perturbed Klein-Gordon equation, Eq. (ff0|), for Sx,k 
then becomes 



(SxkT + 3H(6 X k) 



i 2 2 i o \ 2 

— + m+gip + 3A x x 
or A 



4x$ fe - 2 [g(p 2 x + m 2 xX + X x x 3 ] 
- ZgLpx^Vk- 



5Xk 



(67) 



From this it is clear that the effective mass of the x field 
depends on m x , x 2 an d 9V 2 (whether one uses [(<Sx)/c=o] 2 
or {5x 2 } for the background field x is n °t very impor- 
tant at this qualitative level). In the simple preheating 
case in which metric perturbations are ignored, increas- 
ing m x leads to a sharp decrease in the strength of the 
resonance. (Effectively it increases the A parameter of 
the Mathieu equation (A oc k 2 /a 2 + m^/m^ + 2q), lead- 
ing to a rapid decrease in fxu-) Similarly, a positive A x 
term tends to shut off the x resonance because once the 
variance of the x field becomes large, the 3A x x 2 term acts 
like an effective mass term, hence damping the resonance 
|yj . However in our case we have the extra driving term 
proportional to which also depends on m x and A x . 
Hence it is not obvious that adding these two parame- 
ters will definitely damp the resonance. This remains an 
interesting topic for further research. 



D. Escape route 4: Suppression of Super-Hubble 
Initial Conditions 

The largest resonance effects occur when q ^> 1, that 
is, when the coupling g is much larger than (m/M p \) 2 . 
During inflation, when tp is slowly rolling, such a large 
coupling g would lead to a large effective mass for the x 
field, m x ~ gip 2 S> m 2 . It has been drawn to our atten- 
tion that this large effective m x during inflation might 
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suppress the value of the \ field at the beginning of pre- 
heating on k/a <C H scales, relative to modes in the 
k — > oo limit. Such a suppression of the amplitude 
of the Xk modes in the k — » limit might suggest that 
the large-fc modes would have a "head start" at preheat- 
ing, and would therefore go nonlinear before the super- 
Hubble modes would, effectively damping the resonance 
via backreaction before the super-Hubble modes experi- 
enced significant resonant growth. 

This situation, however, is complicated by several fac- 
tors, each of which deserves further attention: 

• The Floquet indices fik — * rapidly with grow- 
ing k. During the first few inflaton oscilla- 
tions, for multi-field models with large q, modes 
grow due to stochastic resonance. In this case, 
the amplification factor goes as || e Mfc — 1 ~ 
2 exp[— Trk 2 / (2a 2 m 2 y/q e s)], which obviously decays 
rapidly with k. This shows that large-fc modes will 
not grow at a higher exponential rate during pre- 
heating, even if they had a large "head start." 

At later times, when the evolution is closer to ordi- 
nary Floquet theory, the above conclusion remains: 
/Xfc will be larger for the lowest resonance band, at 
small-fc, and drops off rapidly with increasing k. 

• Assuming that long-wavelength modes of interest 
first cross the Hubble scale 50 e-folds before the 
end of inflation, these modes would naively be sup- 
pressed by a factor of a~ 3 / 2 ~ e~ 75 ~ 10~ 30 . 
Even a suppression of Sxk(to) initial conditions for 
the small-fc modes of order 10~ 30 compared with 
large-fc modes is not necessarily significant when 
one considers exponential instabilities at preheat- 
ing: /ife/m ~ 0(1) for fc — > modes when q ^> 1, 
so that small-fc modes could still go nonlinear first 
if Mfe~o > 75/ife_ +00 . Considering that /ij, falls off 
exponentially with fc, this constraint can easily be 
satisfied over large regions of parameter space. 

• Not only are the /ifc much smaller for large fc than 
for small fc, but the widths of the resonance bands, 
Afc, are likewise much narrower for large fc. This is 
crucial since backreaction is controlled by fc-space 
integrals, such as S/ ~ f dkk 2 \Sipik\ 2 , for the field 
fluctuations Sipj. The decrease in the resonance- 
band window at large fc is partially compensated 
for by the larger phase-space volume, k 2 dk factor, 
and therefore requires careful quantitative study. 

• This scenario for suppressing the amplitude of 
super-Hubble modes neglects the quantum-to- 
classical transition, which remains very subtle. 
The standard procedure in nonequilibrium quan- 
tum field theory is to define the initial conditions 

Stfikih) = VV 2fi fe(*o) and [dSipi k /dt] to = 

. This prescription 



has been used in every numerical and nearly ev- 
ery analytical study of preheating. || Note that 
these initial conditions for the start of preheating 
suppress large-fc modes and favor small-fc modes. 
By adopting our scale-invariant initial conditions 
for Stpik in our simulations in Sec. VIA, we have 



thereby neglected this additional enhancement of 
super-Hubble modes. 

Because of the many subtleties involved in the ques- 
tion of initial conditions, this obviously warrants further 
study. It is crucial to note, however, that even if super- 
Hubble modes were suppressed by a large factor during 
inflation, some modes are still very likely to go non- 
linear during preheating, at some scale fc: if preheat- 
ing happens to the matter-field modes at any particular 
scale fc, then the metric perturbations must necessar- 
ily experience similar exponential growth; this is fixed 
by the constraint e quat ion, Eq. (|4l[). And, as empha- 
sized above in Sec. VU, since preheating only ends via 
nonlinear mechanisms, the exponential amplification of 
these $fc modes will continue into the nonlinear regime. 
Whether or not the strongest amplification and earliest 
transition to nonlinearity occurs on super-Hubble or sub- 
Hubble scales, it is bound to happen in some region of 
fc-space during preheating. 



IX. HIERARCHICAL CLASSIFICATION OF 
INFLATIONARY MODELS 

Here we examine a variety of inflationary scenarios for 
possible sensitivity to strong preheating. 

Simple chaotic inflation 

The start of reheating in chaotic inflation is charac- 
terized by very large values of ipo w M p i.™§ In simple 
polynomial models V — A0 2 ™, A is forced to be very 
small due to the COBE CMB results, with A w 1CT 12 for 
n = 2, or m w 10~ 6 M p \ for n = 1. Hence we find the 
result for scalar fields coupled to the inflaton, as in Eq. 



gxlO 



10 



(68) 



as 



iy/Sl k (t a )/2 for all fc. Rj2 



^ § The exact moment when oscillations are said to start, 
and hence the value of the initial amplitude of oscillations, 
ipo, is rather subtle. Using the criterion that it is denned 
by the time when the first of the slow-roll parameters e or n 
become unity, gives ipo = %/2M p i for the quadratic potential 
and ipo = 2\/3Afpi for the quartic potential. Using such large 
values makes the resonance much stronger. To obtain a lower 
limit we use ipo w 0.3M p i, as typical in the literature, see e.g. 
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so that for natural couplings g G [10~ 6 ,10 -2 ] we have 
q € [10 4 , 10 8 ]. In other words, large resonance factors are 
generic in the simplest models of chaotic inflation, when 
scalar fields in addition to the inflaton are included. 

A very important point to notice is that the large q 
arise due to the confluence of the largeness of ipo and the 
smallncss of the inflaton mass during reheating. If either 
of these is given up, the resonances can be expected to 
be naturally much weaker. 

Hence, in a chaotic inflation model based on a potential 
which exhibits a large change in the inflaton effective 
mass between the inflationary phase and the reheating 
phase, the q parameter might easily be several orders 
of magnitude smaller due to the strong curvature of the 
potential at reheating. 

Finally, in the case where the potential exhibits self- 
interaction for the inflaton, e.g. V = \<fi 4 /4:, it has been 
strongly argued that the value A ~ 10~ 12 is unnecessarily 
small, arising from a naive identification of classical per- 
turbations (responsible for the CMB anisotropies) with 
quantum fluctuations, and that a value nearer A ~ 10~ 6 
is more accurate If this is true, then for models 

with a massless quartically self-coupled inflaton coupled 
to distinct bosons, the resonance parameter q would be 
reduced by around 10 6 , since the square of effective fre- 
quency of the inflaton's oscillations, to 2 cx Xtpg, appears 
in the denominator for q. Naturally, this would be a very 
large change, which, for moderate couplings g < 10~ 3 , 
would move the system into the mild or weak resonance 
regimes. This change in A, however, would not neces- 
sarily weaken the resonances for chaotic inflation models 
with a massive inflaton coupled to distinct fields. 

Hybrid and multiple inflation 

General hybrid models can exhibit strong preheating 
with large q parameters 1 85 , and to this extent they will 
suffer strong distortion of the metric perturbation spec- 
trum during any oscillatory phase. What would be inter- 
esting would be to study models in which an oscillatory 
phase with large resonance parameters occurs between 
successive phases of inflation. Once the oscillations have 
finished, the amplified spectrum would be smoothed and 
stretched during the second inflationary phase leaving 
a highly non-trivial imprint on the CMB, assuming the 
second phase did not last too long. 

The extra freedom in these models implies that it 
should be possible to construct models with very simi- 
lar features to those of the broken spectral index genre, 
which are able to fit current CMB and large scale 
structure data better than standard inflationary models 
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FIG. 16. A schematic representation of the hierarchy of in- 
flationary models with regard to their strength of preheating. 
At the top, where q = 0, sits warm inflation. At the bottom, 
where typically q 3> 1, sit the simplest chaotic inflation mod- 
els. Most other inflationary models lie in between and exhibit 
the effects discussed in this paper to varying degrees depend- 
ing on precise details of the couplings between the inflaton 
and other fields. 



Supersymmetric (SUSY) inflation 

Consider perhaps the simplest SUSY theory based on 
the superpotential [[73[: 



W — aStpip — /j, 2 S. 



(69) 



This is the most general form consistent with R-parity 
under which the inflaton singlet, S, transforms as S — > 
e %e S, the superpotential transforms as W — > e l6 W, and 
the product Tpip is invariant. The unbroken supersym- 
metric potential corresponding to this superpotential is: 



V = a 2 \S 2 \(\^ + m 2 ) 



\aipip — /j 2 ] 2 + D-terms . 



(70) 



The £)-terms vanish along the Z?-flat direction \ip\ = \tp*\. 
Inflation occurs for S 3> /Vv^ a t the minimum of the 
potential ((</?) = (Tp) — 0), which is not supersymmetric, 
with V = /i 4 . Quantum corrections break the flatness of 
the potential and cause slow-roll but as long as super- 
symmetry is broken softly the induced curvature will be 
small and the inflaton very light, as in chaotic inflation 
above. If this were to remain the situation, these theo- 
ries would have the same problems as the simple chaotic 
inflation models. However, for S < fi/y^a, V has a new 
minimum which is supersymmetric (V — 0) correspond- 
ing to (S) = 0, (ip) = [ij\fa |73]. Reheating is expected 
to complete during oscillations about this new minimum. 
The amplitude of both S and tp oscillations is of order 
/i/v^a, while their masses are of order m| ss 4/i 2 and 
to 2 ~ 8a/x 2 . Hence, assuming no other fields, the reso- 
nance parameters at the final oscillation stage are: 
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0.1, 



(71) 



(72) 



since a ~ 10 -2 to fit observations |7j|. Hence, in this 
most minimal of supersymmetric models, there is little 
resonance in the primary oscillation phase and in the 
inflaton and Higgs superfields ip. In this sense it is closer 
to the single-field models since it contains only one energy 
scale (/i 2 ) and a single dimensionless coupling, a, which 
is forced to be small. 

However, if we couple for example, another scalar field 
X to S and ip, with coupling g 2 3> a 2 , then the resonance 
parameters would of the order g 2 /a 2 ~ g 2 x 10 , so that 
if g were near strong coupling, the strong resonances we 
have discussed here would occur. However, since these 
extra couplings would typically give mass to other fields, 



(e.g. the heavy right-handed neutrinos of the theory |87|) 
and would play a role in determing the baryon asymme- 
try, building a fully consistent model with g 2 ~ 1 appears 
technically very difficult. 

Hidden sector inflation 



In these string-inspired models |70j, the inflaton in 
the hidden sector only has gravitational couplings to the 
fields of the visible sector and hence the q parameter is 
very small, being of order unity, as occurs for gravita- 
tional waves jl9| . The amplification of metric perturba- 
tions is essentially negligible in these models, as it is in 
the quintessential inflation model of Vilenkin and Peebles 



Warm inflation 

At the opposite end of the scale from simple chaotic 
inflation models is warm inflation and its variants [ ^6| , 
in which inflation occurs due to strong overdamping of 
the inflaton due to couplings to other fields. Reheating 
does not occur due to coherent oscillations but rather 
occurs continuously during inflation with natural exit to 
a radiation-dominated FRW universe. Since there are no 
coherent oscillations to speak of, q = 0, and there is no 
resonant amplification of metric perturbations at all. 



X. CONCLUSIONS AND FUTURE ISSUES 

This paper continues the exploration of metric per- 
turbations' evolution during the post-inflationary epoch 
known as preheating. In spirit it is a continuation and 
development of the study in Paper I [il^j , but with sig- 
nificant advances in conceptual and numerical analysis of 



the problem. We have identified a number of important 
issues that distinguish the evolution of metric perturba- 
tions during the violent preheating epoch (q 3> 1) from 
metric evolution in old theories of reheating. Primary 
among these is the resonant amplification of both super- 
and sub- horizon modes of <&, the gauge-invariant scalar 
metric perturbation. 

We demonstrate explicitly that such amplification does 
not violate causality - it is merely a reflection of the 
perfection of the coherence of the inflaton condensate at 
the very end of inflation, which is correlated on vastly 
super-Hubble scales in standard models of inflation. 

Such resonances depend heavily upon the presence of 
multiple matter-fields coupled via non-gravitational in- 
teractions. In the multi-field case, entropy perturbations 
are generic, and although they may be small during in- 
flation, they are resonantly amplified during strong pre- 
heating. This leads to the destruction of the conservation 
of traditional quantities such as the Bardeen parameter 
C, traditionally used to transfer the metric perturbation 
4> from inflation to photon decoupling. This implies that 
direct numerical integration of the governing equations 
must be used. 

The wide resonances at g > 1 destroy the scale- 
invariance of the inflationary power spectrum and could 
hugely over-produce large-angle anisotropies in the Cos- 
mic Microwave Background if left undamped after pre- 
heating. The resonances can be so strong as to force 
the perturbations to go nonlinear, requiring study of the 
full, unlinearized Einstein field equations (see (25| for an 
exciting start in this direction in the single- field case). 
One issue that is likely to be very important in under- 
standing the full extent of amplification of the super- 
Hubble modes is the break-up of the inflaton condensate. 
As second-order effects become important, we expect on 
general grounds that the perfect coherence of the oscilla- 
tions of the inflaton condensate will be lost, and the co- 
herence length of the homogeneous mode will be reduced. 
This will tend to shut off the super-Hubble resonances, 
following essentially from the induced super-Hubble in- 
homogeneity. 

We wish to emphasize that this powerful growth in 
the metric perturbations is robust and cannot be swept 
away or removed by fine-tuning the initial value of <f> 
during inflation. The primary reason why fine-tuning 
fails is that the resonances of preheating end typically 
only through backreaction effects, which require at least 
study of second-order metric perturbation effects and a 
thorough understanding of mode-mode coupling. This is 
the first time inflationary models have been faced with 
this complication, bringing them nearer in spirit to the 
nonlinear field theories used in topological defect stud- 
ies. This mode-mode coupling tends to re-establish the 
scale-invariance of the power spectrum, but this time 
with a spectral index which is determined by the non- 
linear structure of the Einstein field equations, in anal- 
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ogy with the situation with turbulence and the Navier- 
Stokes equations. Memory of the initial spectral index 
is lost and the perturbations tend to an active, incoher- 
ent state, again mimicking topological defects. This may 
have profound implications for the small-scale CMB and 
in particular for the secondary Doppler peaks, as a means 
for distinguishing between inflation and defect theories, 
assuming the mode-mode coupling is preserved on rele- 
vant scales until decoupling. 

As a means for quantifying the amplification of per- 
turbations, we have introduced here the concept of the 
time to nonlinear ity, t n \, which provides a robust (though 
non-unique) and useful generalization of the Floquet in- 
dex Hk for expanding universes. It quantifies the time at 
which the linearized equations break down. 

The metric preheating we have developed here, despite 
being based on a simple model, yields very strong de- 
viations from the old theory of perturbative evolution. 
These deviations depend crucially on a number of con- 
ceptual points, highlighted in Sec. Ill, related to the ex- 



istence of multiple, non-gravitationally interacting fields 
at preheating. While certain regions of preheating q- 
parameter space appear to be observationally ruled out 
by the CMB, there appear a number of interesting and 
even exciting mechanisms for making the strong reso- 
nances acceptable. Secondary phases of inflation are ex- 
pected to damp the amplified modes and stretch them. 
This allows one to put more freedom into the original 
spectrum while bringing the overall amplitude down to 
an observationally acceptable level. 

Alternatively, the resonance may be significantly 
weaker if the fields into which the inflaton is decaying 
have strong self-interaction, or are fermionic. The decay 
into fermions is particularly interesting since it is con- 
strained by the Pauli exclusion principle, which limits 
the amplification of $. While we have discussed the ba- 
sic mechanisms at work, both of these aspects deserve 
closer attention. 

For the theory of preheating itself, one of the impli- 
cations of studying the complete system of equations in- 
cluding metric perturbations is the appearance of rescat- 
tcring - the production of quanta due to couplings be- 
tween the fluctuations of different fields - even at linear 
order. Previously this had been thought Q to be a purely 
nonlinear process. On smaller scales, the amplification of 
scalar metric perturbations strongly enhances primordial 
black hole formation and gravitational wave emission via 
rescattering. 

We note that all inflationary models are definitely not 
born equal with regard to these resonances. Simple 
chaotic inflation models, with their large amplitude oscil- 
lations and small effective mass, yield extremely strong 
resonances. Alternatively, models with big changes to 
the effective mass and/or small amplitude oscillations or 
weak couplings between the inflaton and other fields, ex- 
hibit much weaker resonances. 



To conclude, we believe it is important to analyze the 
limitations of the work presented here. Since preheating 
is a problem involving multi-field, non-equilibrium, semi- 
classical quantum gravity on a dynamic background, it 
is clear that the essentially classical model we have used 
is not ideal. In addition, we have uncovered a serious 
ambiguity related to the definition of the homogeneous 
components of the non-inflaton fields. The standard def- 
inition in the metric perturbation literature ignores the 
contribution of <&(k = 0) to the evolution of homoge- 
nous background quantities. Since the energy in this 
mode grows very rapidly, we have argued that its in- 
fluence on the background should be included as a first 
step to studying the full backreaction problem. 

Other aspects of preheating that should be addressed 
in future work include: 

(i) Can one formulate a consistent quantum and 
non-equilibrium framework which includes the met- 
ric perturbations? In particular, does the Schwinger- 
Keldysh closed-time-path (CTP) formalism @,|^j89 
carry through to this more complex setting? Further, 
can one reconcile the standard relativistic definitions of 
background quantities with the order-parameter | p4) and 
other field-theoretic definitions? 

(ii) Backreaction issues: can one go beyond the simple 
linearized equations in a realistic, analytical manner? Is 
there a relatively simple and consistent Hartree or large- 
N approximation for Eqs. (H^-[ll])? 

(iii) The transfer function to decoupling and thermal- 
ization: how does thermalization proceed and how do 
metric perturbations evolve through this phase? How 
does the viscosity arising from finite temperature effects 
alter perturbation evolution? |4^j59| , |90| , |9l[ | 

(iv) Quantum gravity corrections: stochastic backreac- 
tion effects arc expected to be strong during preheating. 
This implies a stochastic evolution of background quan- 
tities such as the scale factor and zero modes |39|. The 
evolution of the metric is determined by a Langevin-type 
equation which describes the non-equilibrium dynamics 
of the gravitational field. The effects of these stochastic 
fluctuations have already been analyzed phenomenologi- 
cally in the context of preheating without metric pertur- 
bations [34|j3^j3l[| , and found generically to enhance par- 
ticle production. We have similar expectations on very 
general grounds for the case in which metric fluctuations 
are included. (See also p2| .) Almost surely going to the 
(white-noise) stochastic limit will not help reduce the am- 
plification of metric perturbations, though it will change 
the wavelength-dependence of the amplification. 
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